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FOREWORD 


The  current  work  is  in  response  to  the  demand  for  a  computational  tool  that  can  be  used  for 
rapid  and  inexpensive  aerodynamic  analysis  of  missile-type  geometries.  Although  it  is  possible  to 
evaluate  missile  geometries  using  a  full  Navier-Stokes  computational  fluid  dynamics  procedure 
(CFD),  both  budgetary  and  time  constraints  often  preclude  this  type  of  in-depth  analysis,  "pie 
Zeus’^  Flowfield  Analysis  Took  was  developed  in  order  to  minimize  the  cost  and  turnaround  time 
associated  with  aerodynamic  analysis  of  missile  geometries.  The  tool  can  be  used  to  either  perform 
a  flowfield  analysis,  or  alternatively,  to  generate  computational  domains  for  a  full  Navier-Stokes 
solver. 


Initial  funding  for  the  work  was  provided  by  the  Modeling  and  Simulation  Task  of  the  Air 
and  Surface  Weaponry  Technology  Program  funded  and  managed  by  the  Office  of  Naval 
Research,  Code  351.  Additional  funding  was  provided  by  George  Long  through  the  Standard 
Missile  Program,  PMS422.  The  author  would  like  to  express  his  appreciation  to  both  Gil  Graff 
and  George  Long  for  their  support  during  this  project. 

This  report  has  been  reviewed  by  Mr.  Thomas  A.  DelGuidice,  Head,  Weapons  Integration 
and  Technology  Branch;  and  Mr.  Bob  Steigler,  Head,  Combat  Systems  Safety  and  Engineering 
Division. 
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CHAPTER  1:  INTRODUCTION 


1.1  OVERVIEW  OF  ZEUS^ 

In  the  past  decade,  significant  advances  have  been  made  in  the  field  of  computational 
fluid  dynamics  (CFD).  These  include  improvements  in  computational  speed,  numeric^ 
algorithm  efficiency,  and  modeling  (turbulence  and  combustion).  Because  of  this 
progress,  it  is  now  commonplace  to  see  accurate  CFD  solutions  to  very  complex,  large 
scale,  three-dimensional  flowfields.  It  should  be  noted,  however,  that  a  number  of 
limitations  remain.  First,  existing  CFD  solvers  are  quite  complex  and  extensive  training  is 
required  prior  to  generating  a  solution.  The  user  must  be  trained  not  only  in  basic 
aerodynamics,  but  also  in  computational  domain  development  and  numerical  methods. 
Second,  the  conqiutational  costs  are  still  extensive.  Complex  cases  can  still  take  days, 
weeks,  or  even  months  of  expensive  workstation  or  supercomputer  time.  Because  of  these 
limitations,  parametric  type  design  studies  are  often  infeasible,  and  therefore,  the  use  of 
CFD  as  an  engineering  design/development  tool  has  been  limited. 

The  goal  of  the  current  work  is  to  significantly  reduce  both  the  costs  and  turnaround 
time  associated  with  CFD  solutions  of  missile  t^  geometries  with  arbitrary  control 
surfaces.  The  current  work  began  with  the  ZEUS  flowfield  solver^  '’  ’  ®.  The  ZEUS  code 
is  a  supersonic  space  marching  Euler  solver  capable  of  computing  both  internal  and  external 
flowfields.  The  current  work  involved  modifying  the  ZEUS  code  to  make  it  more 
amenable  to  solving  for  the  forces/moments  on  missile  type  geometries  with  arbitrary 
control  surfaces.  A  number  of  limitations  in  the  original  code  motivated  the  current  work:  ■ 

1)  For  each  geometry  under  consideration,  the  user  must  develop  a  Fortran  source 
code  that  provides  a  geometric  description  of  the  missile  surface.  In  addition, 
code  that  specifies  die  control  surfaces  and  switching  logic  for  the  boundary 
conditions  must  be  generated. 

2)  The  user  must  be  familiar  with  a  number  of  separate  codes  in  order  to 
implement  a  ZEUS  run. 

These  include  a  separate  code  for  each  of  the  following  tasks: 

a)  Cone/Blunt  -  Generate  an  initial  profile  by  running  one  of  the  two 
codes. 

b)  Convert  -  Convert  the  initial  profile  to  a  form  readable  by  ZEUS 

c)  ZEUS  -  Solves  the  Euler  form  of  the  equations  of  motion. 

d)  BL  -  Performs  a  boundary  layer  analysis  of  the  inviscid  flowfield. 

Since  the  user  must  understand  all  of  the  inputs/outputs  from  each  of  the  above 
codes,  this  significantly  increases  the  amount  of  training  required. 
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3)  The  setup  time  required  before  executing  the  code  can  be  significant  (typically 
on  the  order  of  hours/days).  In  addition,  if  parametric  studies  are  desired,  a 
separate  set  of  input  files  must  be  generated  for  each  case  under  consideration. 
This  is  both  a  time  consuming  and  error  prone  process  when  a  large  number  of 
cases  are  desired. 

4)  Inherent  limitations  (space  marching)  in  the  ZEUS  code  prevent  the  solution  of 
any  flowfield  with  subsonic  regions. 

The  modifications  made  to  the  ZEUS  code  focused  on  the  following  four  primary 
areas  of  improvement: 

1)  Minimize  the  required  training  -  This  was  done  by  developing  a  graphical  user 
interface,  modifying  the  grid  generation  capabilities,  and  developing  simple  to 
use  pre-  and  post-processing  features.  The  GUI  allows  for  most  of  the 
complexity  of  the  code  to  be  hidden  from  the  user.  Instead  of  developing 
Fortran  modules  to  specify  the  missile  and  control  surface  geometries,  die 
current  tool  (Zeus‘S  uses  a  series  of  point-and-click  GUI  menus  to  setup  the 
required  geometry  files.  All  of  the  complexities  of  generating  a  three- 
dimensional  grid  and  blending  control  surfaces  to  the  missile  have  been  hidden 
from  the  user.  The  GUI  interface  also  significandy  simplifies  the  post¬ 
processing  of  the  results. 

2)  Reduce  the  setup/execution  time  -  The  interface  reduces  the  setup  time  required 
to  a  few  minutes  vs.  many  hours/days  for  the  original  ZEUS  code.  In  addition, 
the  user  can  setup  the  code  to  sweep  through  a  variety  of  parameters  (Mach 
number,  angle  of  attack,  angle  of  yaw,  and  roll  angle)  by  simply  specifying  an 
initial,  final,  and  incremental  value  for  each  of  the  parameters.  All  of  the 
complexity  involved  with  generating  numerous  input  files  has  been  hidden  from 
the  user.  This  feature  allows  for  parametric  ty^  studies  with  no  additional 
setup  time  required. 

3)  Minimize  die  costs  -  The  costs  associated  with  computational  fluid  dynamic 
solutions  arc  often  prohibitively  high,  and  therefore,  prevent  the  application  of 
CFD  during  the  design/analysis  phase  of  missile  development.  The  two 
primary  influences  on  these  costs  are  personal  and  computational  requirements. 
The  personal  costs  are  inflated  because  you  must  have  someone  with  extensive 
training  performing  a  very  time  consuming  task.  Zeus"^  reduces  the  personnel 
costs  by  minimizing  the  required  training  and  reducing  the  time  required  to 
setup  a  run.  The  second  factor,  computational  cost,  is  high  because  the  ZEUS 
code  was  developed  for  use  on  expensive  Unix  workstations.  These  costs  were 
reduced  by  porting  the  ZEUS  code  to  less  expensive  personal  computers  and  by 
developing  a  GUI  that  runs  imder  a  widely  available  operating  system 
(Windows  NT  and  Windows  95). 

4)  Extend  the  range  of  applications  for  the  code  -  There  are  a  number  of  cases 
where  either  the  ZEUS  code  caimot  be  applied,  or  where  the  results  may  be 
suspect  due  to  the  approximate  nature  of  the  solution  procedure  (space  marching 
Euler).  For  these  cases,  the  Zeus^  tool  includes  a  grid  generator  which  will 
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output  a  three-dimensional  grid  and  the  appropriate  boundary  conditions  which 
can  then  be  applied  as  input  for  the  GASP*  full  Navier-Stokes  flowfield  solver. 
It  should  be  noted  that  generation  of  complex  three-dimensional  missile  grids 
with  control  surfaces  often  took  days  to  develop.  The  Zeus"^  tool  reduces  this 
to  a  matter  of  minutes  on  a  personal  computer. 


1.2  DCX:UMENTATION 

The  documentation  for  the  Zeus^  tool  is  presented  in  one  of  two  forms:  online  help  or 
a  printed  user’s  manual.  It  should  be  noted  that  the  data  contained  within  the  two  forms  are 
identical. 

To  access  the  online  help  press  the  'F 1 '  key  during  execution  of  the  Zeus**  code.  If 
the  user  does  not  understand  one  of  the  requested  inputs,  pressing  the  'FI'  key  will  display 
the  help  file  for  that  specific  topic.  Alternatively,  the  user  can  search  through  the  help 
documentation  online  by  clicking  on  'Help'  from  the  main  view  screen. 

A  hard  copy  of  the  Zeus**  user’s  manual  can  be  obtained  by  executing  the  code 
(double  clicking  on  ZeusPP.exe)  and  selecting  “Print  User’s  Manual”.  Alternatively,  tte 
mannal  may  be  printed  by  selecting  “Help”  and  then  “Print  User’s  Manual”  from  the  main 

window. 

Note:  The  manual  is  provided  in  Adobe  pdf  format.  The  Adobe  Acrobat  reader 
may  be  obtained  from  Adobe’s  home  page  at  “http://www.adobe.com/acrobat”. 
Alternative  forms  of  the  User’s  manual  can  be  supplied  upon  request  if  the  pdf 
format  cannot  be  viewed/printed. 


1.3  SYSTEM  REQUIREMENTS 


Processor: 
Operating  System: 

Plotting  Package: 
Hard  Disk  Space: 


Intel  486  compatible  (233  MHz  Pentium  II  or  higher  recommended) 
Microsoft  Windows  NT  4.0  (recommended),  Windows  95,  or 
Windows  98 
Tecplot  v7.0  (Optional) 

Approximately  3  Megabytes  for  the  Zeus**  code  and  associated  files. 
The  size  of  the  output  data  generated  by  Zeus**  is  highly  case 
dependent  (grid  resolution,  #  of  sweeps,  etc. ..)  and  can  range  from  a 
few  megabytes  to  a  few  hundred  megabytes. 


1.4  CONTACT  THE  AUTHOR 

E-mail  requests  for  a  copy  of  the  code  to  the  author  at  the  following: 

Dr.  David  F.  Robinson  (Code  G72) 

Naval  Surface  Warfare  Center 
17320  Dahlgren  Rd. 

Dahlgren,  VA  22448-5100 
robinsondf@nswc.navy.mil 


3 


NSWCDD/TR-98/147 


Voice:  (540)653-3954 
Fax:  (540)  653-8268 


1.5  EXECUTE  THE  CODE 

Execution  of  the  code  begins  by  double  clicking  on  the  program  'ZeusPP.exe’, 
entering  a  name  for  the  current  case,  and  then  clicking  on  ’start  New  Case'. 

Note:  1)  The  name  entered  cannot  contain  blank  spaces. 

2)  A  tutorial  for  setting  up  and  executing  a  new  case  is  provided  in 
Appendix  A 

3)  The  code  must  be  executed  from  the  main  ZeusPP  directory.  To  execute 
one  of  the  sample  cases  provided,  copy  the  file  from  the  ’samples'  directory 
to  the  ZeusPP  main  directory. 

To  open  an  existing  case  (*.zpp)  click  on  'ZeusPP.exe'  and  select  'Open  Existing 
Case'.  From  the  main  menu  click  on  'File/Open'  and  select  the  desired  case.  Alternatively, 
simply  double  click  on  the  file  ('filename.zpp'). 
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CHAPTER  2:  MAIN  WINDOW  CONTROLS 


2.1  GEOMETRY  SKETCH 

As  the  geometry  of  the 
dimensional  sketches  of  the  side 
There  are  a 
upper  left 

below. 


. ^ —  - . - . . . - . 

FIGURE  1.  MAIN  VIEW  SCREEN 


(see  example  below), 
which  is  listed  in  the 


missile  and  its  control  surfaces  are  specified,  two- 
and  frontal  views  are  displayed  in  the  main  view  screen 
I  variety  of  methods  for  manipulating  the  sketch,  each  of 
comer  of  the  main  view  screen  and  discussed  in  detail 


•  Ctrl-R  -  Pressing  the  'Control'  and  'R'  keys  simultaneously  will  regenerate  the 
geometry  display  shown  in  the  main  window.  This  should  be  done  any  time  the 
missile  geometry  is  modified  and  an  updated  view  is  desired.  Note  that  this  view 
is  for  display  purposes  only,  and  therefore,  to  limit  computational  effort,  a  low- 
resolution  representation  of  the  surface  geometry  is  provided.  This  is  not 
representative  of  the  computational  domain  employed  during  the  solution 
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procedure.  If  desired,  the  resolution  of  the  display  can  be  increased/decreased  by 
modifying  the  number  of  points  in  the  View/Options  section. 

•  +/. .  Pressing  the  '+'  and  keys  steps  the  views  circumferentially  through  each  of 

the  zone  edges.  '  +  ’  shows  die  next  and  shows  the  previous  circumferential 
view. 

•  n/p-  Pressing  the  'n'  and  'p'  keys  steps  the  front  view  axially  down  the  missile.  '  n ' 

and  ’p'  go  to  the  next  and  previous  axial  location,  respectively.  Pressing  'N'  or 
'P'  has  the  same  effect,  however,  larger  axial  steps  are  taken.  The  current  axial 
location  is  shown  in  the  lower  left  comer  of  the  main  view  screen  (z=0.275). 

•  A/a-  Pressing  the  'A'  or  'a'  key  redraws  all  of  the  axial  stations  and  circumferential 

planes. 

•  Ctrl-B  -  Pressing  the  'Control'  and  'B'  keys  simultaneously  will  display  the  axial 

stations  where  the  inviscid  grid  was  output.  These  stations  are  shown  as  blue  tic 
marks  along  the  centerline  of  the  missile  geometry.  This  feature  cannot  be 
employed  until  after  an  inviscid  run  has  been  completed. 

Note:  The  viscous  boundary  layer  solver  requires  the  inviscid  Euler  solution  as  a 
boundary  condition.  Therefore,  the  inviscid  solution  must  be  performed 
before  running  tiie  viscous  code.  To  obtain  an  accurate  viscous 
approximation  an  adequate  number  of  inviscid  profiles  must  be  provided 
such  that  all  relevant  geometric  features  are  spati^y  resolved.  The  inviscid 
profiles  are  written  out  at  the  user-specified  intervals  provided  in  the 
boundary  layer  inputs  section. 

•  'z'/'Z'  -  Pressing  the  'z'  or  'Z'  keys  zooms  in/out  on  die  main  view  screen.  The 

amount  of  scaling  applied  to  the  original  sketch  is  shown  in  the  upper  right  hand 
comer.  To  return  to  &e  original  settings,  press  Ctrl-R  to  regenerate  the  sketch. 

•  Arrow  keys  -  The  arrow  keys  are  used  to  translate  the  sketches  in  the  main  window. 


2.2  MENU  OPTIONS 

File 

View 

Geometry  Setup 

Inputs 

Execution 

Post-Processor 

Window 

Help 


Controls  file  I/O  and  printing. 

Toolbar  views  and  runtime  options. 

Define  the  missile  and  control  surface  geometries  as  well  as 
the  computational  domain. 

Input  free-stream  conditions,  runtime  parameters,  reference 
conditions,  ...etc. 

Executing  the  code 
Post-Processing 

Layout  and  orientation  of  windows 

Help  controls  and  'About'  information  (author,  version  #, 
and  contact  information) 
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2.3  MISCELLANEOUS  OPTIONS 


nGURE2.  PROGRAM  CONTROL  PARAMETERS 


•  Number  of  Points  -  The  sketch  shown  in  the  main  window  is  a  low-resolution 

representation  of  the  current  missile  geometry.  To  increase/decrease  the  resolution 
of  the  geometric  display  simply  increase/decrease  the  appropriate  points  parameter 
(axial,  radial,  or  circumferential).  Increasing  the  resolution  will  increase  both  the 
time  to  display  the  sketch  as  well  as  the  size  of  the  output  data  file  (*.zpp). 

•  Save  sketch  -  The  output  data  file  for  the  Zeus"^  code.(*.zpp)  contains  all  of  the 

input  parameters  as  well  as  the  coordinates  for  the  geometry  sketch  shown  in  the 
main  window.  To  minimize  the  size  of  the  output  files  (*.zpp)  de-select  the  save 
sketch  on  exit'  dialog  box  shown  above.  This  will  significantly  decrease  the  size 
of  the  output  file  (approximately  3k  vs.  1-2  meg.);  however,  when  the  case  is 
restarted,  the  geometry  sketch  must  be  regenerated.  If  the  size  of  the  output  file  is 
not  important  then  it  is  recommend  that  this  option  be  selected  to  minimize  the 
computational  effort  required  at  startup. 

•  Tecplot  Parameters  -  A  number  of  Tecplot  macro  files  (contour  plots,  movie  files, 

. .  .etc.)  are  generated  by  the  Zeus‘S  code  in  order  to  simplify  post-processing.  The 
following  parameters,  used  in  the  macro  files,  are  specified  in  the  View/Options 
section. 

X  • ,  X,  -  The  minimum/maximum  x-axis  range  for  the  movie  files. 

Y,^n»  -  The  minimum/maximum  y-axis  range  for  the  movie  files. 

Flood/Lhie  -  Specify  either  flood  or  line  drawing  for  the  contour  plots. 

#  of  contours  -  Specify  the  number  of  contour  lines  to  be  drawn. 

•  Pause  -  Zeus"'^  runs  a  number  of  external  programs  from  the  command  prompt.  If 

the  'Pause  after  each  run'  feature  is  selected  then  the  command  prompt  will  remain 
active  (i.e.,  paused)  and  prompt  the  user  to  'press  any  key'  before  exiting.  If  the 
pause  feature  is  false,  the  Zeus"^  code  will  automatically  close  the  command 
prompt  upon  termination  of  the  code. 
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CHAPTER  3;  GEOMETRY  SETUP/EVALUATION 


MENU  OPTIONS 


Grid  Parameters 


Missile  Sections 
Control  Surfaces 
Grid  Generation 


Specify  the  grid  parameters  for  both  the  Euler  and  Navier-Stokes 
grids.  (#  of  grid  points,  clustering,  #  of  zones,  and  location  of  outer 
boundary). 

Define  the  missile  geometry 

Define  the  control  surface  geometry 

Generate  the  Euler  and  Navier-Stokes  computational  domains. 


3. 1  COMPUTATIONAL  DOMAIN 

3.1.1  Grid  Parameters 

Recall  the  Zeus‘S  tool  has  two  primary  functions.  The  code  serves  not  only  as  an 
approximate  flow  solver,  but  also  as  a  Navier-Stokes  grid  generator.  Therefore,  two 
separate  sets  of  grid  parameters  are  specified  when  setting  up  the  computational  domain.  In 
the  first,  the  user  inputs  the  parameters  for  the  computational  domain  that  will  be  used 
during  the  Euler  and  boundary  layer  solutions  of  the  flowfield.  In  the  second,  the  input 
parameters  are  used  to  generate  the  Navier-Stokes  grid  and  appropriate  boundary 
conditions. 

Mesh  Dimensions 

•  Number  of  zones  -  The  number  of  zones  is  equal  to  the  maximum  number  of 
control  surfaces  at  a  given  axial  location  (i.e.,  if  you  have  4  canards,  3  dorsal 
fins,  and  5  tails  fins  then  you  must  enter  '5'  as  the  number  of  zones). 

•  No  control  surfaces  -  Select  this  box  if  you  are  considering  a  case  with  no 
control  surfaces. 

•  Number  of  radial  cells  -  The  number  of  cells  located  radially  between  the  body 
surface  and  the  outer  boundary  of  the  computational  domain. 

•  Default  number  of  circumferential  cells/zone  -  Changing  this  parameter  will 
change  the  number  of  circumferential  cells/zone  for  all  zones. 

•  Number  of  cells/zone  -  The  number  of  cells/zone  refers  to  die  circumferential 
direction  (i.e.,  between  zone  edges)  and  does  not  have  to  be  the  same  for  all 
zones.  Change  an  individual  zone  by  clicking  on  the  desired  zone  and  entering 
a  new  value. 
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Mesh  Clustering 

The  flowfield  solver  typically  does  not  require  mesh  clustering  because  there  is  no 
boundary  layer  to  resolve.  However,  in  certain  cases,  mesh  clustering  (radial  and 
circumferential)  can  be  used  to  more  accurately  capture  specific  flowfield  features 
(vortices,  shocks,  etc).  Note  that  mesh  clustering  is  required  for  the  Navier-Stokes 
grid  in  order  to  resolve  the  boundary  layer. 

•  Radial  clustering  -  Specify  the  grid  clustering  between  the  body  surface  and  the 
outer  boundary.  A  value  of  1 .0  provides  equal  spacing  and  a  value  of  0.0 
provides  maximum  clustering  at  the  surface. 

Note:  1)  Two  sets  of  radial  clustering  are  required.  These  correspond  to  the 
radial  clustering  at  the  initial  and  final  (nose  and  tail)  axial  stations.  This 
provides  the  ability  to  have  more  radial  clustering  at  the  nose  than  at  the 
tail  of  the  missile,  or  vice-versa.  The  radial  clustering  parameter  is 
interpolated  linearly  between  the  two  values  as  you  progress  axially 
down  the  missile. 

2)  Clicking  the  'Boundary  Layer  Thickness'  button  provides  an  estimate 
of  the  boundary  thickness  (laminar  and  turbulent)  at  a  given  axial 
station.  The  values  are  calculated  using  flat  plate  boundary  layer 
assumptions  and  can  be  used  to  estimate  the  amount  of  radial  clustering 
rpqiiirp.fl  for  a  viscous  grid.  This  calculation  requires  a  umt  Reynolds 
number  and  an  axial  location.  Initial  spacing  for  the  viscous  grid  is 
typically  taken  as  l/30th  of  the  estimated  boundary  layer  thickness. 

•  Circumferential  clustering  -  Specify  the  clustering  at  the  zone  edges.  A  value  of 
1.0  provides  equal  spacing  and  a  value  of  0.0  provides  maximum  clustering  at 
the  zone  edges  (i.e.,  circvunferentially  between  control  surfaces). 

Outer  Boundary 

•  Distance  to  outer  boundary  -  If  the  outer  boundary  of  the  computational  domain 
is  not  tracking  a  shock  (i.e.,  Free-Stream  starting  condition)  then  the  radial 
distance  from  the  body  surface  to  the  outer  boundary  is  specified  in  terms  of 
body  radii  (typically  10-20  body  radii  for  an  inviscid  solution). 

•  Growth  of  outer  boundary  -  In  certain  cases  it  is  desirable  for  the  radial  distance 
between  the  body  surface  and  the  outer  boundary  to  increase  as  you  march 
axially  down  the  missile.  This  parameter  specifies  the  angle  (in  degrees)  which 
the  outer  boundary  grows  relative  to  the  body  surface. 

Note:  The  above  parameters  are  used  during  the  Navier-Stokes  grid 
generation.  They  are  also  used  during  the  Euler  solution  if  the  outer 
boundary  does  not  track  the  conical/bow  shock.  The  type  of  outer 
boundary  applied  by  the  flowfield  solver  is  specified  in  the  Integration 
Parameters  section. 
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Copy  Parameters 

•  Select  'Copy  Parameters'  button  to  copy  the  grid  parameters  from  the  alternate 
grid  (i.e.,from  the  Navier-Stokes  grid  if  you  are  in  the  Eider  grid  setup,  or 
vice-versa) 


3.2  MISSILE  GEOMETRY 
3.2. 1  Select  Geometry  Type 

The  first  step  in  setting  up  the  geometry  is  to  define  all  of  the  sections  that  comprise 
the  missile  surface  (i.e.,  not  control  surfaces).  First,  select  the  desired  type  of  nose  and 
then  click  on  'Add  New  Section'.  You  will  be  prompted  for  the  geometric  parameters  used 
to  generate  the  desired  type  of  nose  section.  After  the  nose  section  has  been  specified, 
repeat  the  above  process  for  any  remaining  sections  (i.e.,  afterbody,  boat-tail,  etc. . .).  Note 
that  the  geometry  sections  must  be  entered  in  the  order  in  which  they  occur  on  the 
geometry. 


HGURE  3.  MISSELE  GEOMETRY  SELECTION 


It  should  be  noted  from  the  above  figure,  that  there  are  a  variety  of  possible  missile 
sections.  A  sample  sketch,  a  description  of  the  required  parameters,  and  the  generating 
equations  are  shown  below  for  each  of  the  available  geometry  sections.  In  all  of  the 
geometric  formulas  discussed  below,  'z'  is  the  axial  distance  from  the  start  of  the  given 
section. 
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3.2.2  Conic 


nGURE4.  CONIC  SECTION  PARAMETERS 


A  variety  of  geometries  can  be  described  using  die  conic  input  section.  These 
include  those  shown  in  the  sketch  above,  namely,  a  nose  cone,  a  cylindrical  section,  and  a 
boat-tail  flare.  The  parameters  required  to  generate  a  conic  section  are  shown  m  the  above 
figure. 

3.2.3  Tangent  Ogive 


nOURES.  TANGENT  OGIVE  PARAMETERS 

A  sample  sketch  of  a  tangent  ogive  is  shown  above  along  with  the  parameters 
required  to  generate  the  section.  Note  from  the  sketch  that  the  missile  section  is  tangent  to 
the  horizontal  at  the  aft  axial  station,  hence  the  name,  tangent  ogive.  The  geometry  is 
generated  using  the  following  formula’®  along  with  the  user  specified  parameters  R  and  L, 
where  R  is  the  base  radius  and  L  is  the  length  of  the  section. 


NSWCDD/TR-98/147 

r  =  ^Jp-  +{z-Ly  +{R-p) 
R-  + 1} 


3.2.4  Secant  Ogive 


HGURE  6.  SECANT  OGIVE  PARAMETERS 

A  sample  sketch  of  a  secant  ogive  is  shown  above  along  with  the  parameters 
required  to  generate  the  section.  The  geometry  is  generated  using  the  following  formula'® 
along  with  the  user  specified  parameters  L,  R,  and  p ,  where  L  is  the  length  of  the  section, 
R  is  the  base  radius,  and  p  is  the  radius  of  curvature. 


r  =  -)j  p~  +  {z-h)'  +k 
h  =  pcosa 
k  =  psina 


a  =  y-d 

0  =  tan”' 

7  =  cos"' 

> 
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FIGURE?.  VON  KARMAN  OGIVE  PARAMETERS 


A  sample  sketch  of  a  Von  Karman  ogive  is  shown  above  along  with  the  parametere 
required  to  generate  the  section.  The  geometry  is  generated  using  the  following  fonnula 
along  with  the  user  specified  parameters  R  and  L,  where  R  is  the  base  radius  and  L  is  the 
length  of  the  section. 


0  =  cos  ^(l-2z) 
C  =  0 


3.2.6  Haack  Series 


FIGURE  8.  HAACK  SERIES  NOSE  PARAMETERS 
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A  sample  sketch  of  a  Haack  series  is  shown  above  along  with  the  required  input 
parameters.  The  geometry  is  generated  using  the  following  formula’®  along  with  the  user 
specified  parameters  R  and  L,  where  R  is  the  base  radius  and  L  is  the  length  of  the  section. 


0  =cos  ’(l-2z) 


3.2.7  Power  Series 


HGURE  9.  POWER  SERIES  NOSE  PARAMETERS 


A  sample  sketch  of  a  Power  series  section  is  shown  hbove  along  with  the  required 
input  parameters.  The  geometry  is  generated  using  the  following  formula’®  along  with  the 
user-specified  parameters  R,  L,  and  n,  where  R  is  the  base  radius,  L  is  the  length  of  the 
section,  and  n  is  the  power  series  exponent. 

r=z" 
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3.2.8  Spherical  Nose 


HGURE  10.  SPHERICAL  NOSE  PARAMETERS 

A  sample  sketch  of  a  spherical  nose  section  is  shown  above  along  with  the 
parameters  required  to  generate  the  section.  For  the  spherical  nose,  there  will  be  a  small 
region  where  the  flow  is  subsonic.  Because  the  Euler  solver  employs  a  marching 
procedure,  the  Zeus  code  cannot  be  apphed  in  this  subsonic  region.  Therefore,  the  blunt 
body  solver  must  be  apphed  over  the  entire  subsonic  region.  The  outflow  plane  froin  the 
blunt  body  solver  can  then  be  taken  as  the  inflow  plane  for  the  Euler  solver.  The  input 
parameters  for  the  blunt  body  solver  are  specified  in  the  Integration  Parameters  section. 

3.2.9  r.ro.ss-Sectional  Geometry 

For  each  missile  section  defined,  a  cross-sectional  geometry  must  be  specified.  If  a 
cross-section  is  not  designated  the  default  value  of  'circular'  will  be  used. 
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•  Cross-Section  -  Select  the  desired  cross-sectional  geometry. 

•  Comer  Bluntness  -  Specify  a  value  between  zero  (no  blunting)  and  0.999 
(completely  blunted  comers)  for  the  Rectangular  cross-sectional  geometry. 

•  Blending  -  It  is  often  desirable  to  have  varying  cross-sectional  geometries  for 
different  sections  of  the  missile  (i.e.,  circular  nose  cone  and  elliptical 
afterbody).  For  these  cases,  the  cross-sectional  geometry  caimot  be  abraptly 
changed,  but  rather,  must  be  blended  axially  between  the  two  sections.  If 
blending  is  desired,  then  select  the  'Blend  w/  previous  section'  box  and  specify 
a  blending  distance  in  ft/m. 


3.3  CONTROL  SURFACE  GEOMETRY 
3.3.1  Select  Geometry  Type 

After  the  missile  geometry  has  been  defined,  the  next  step  is  to  add  the  control 
surfaces.  If  no  control  surfaces  will  be  examined,  then  the  'No  Control  Surfaces' 
parameter  must  be  selected  in  the  grid  parameter  section.  Select  the  desired  type  of  control 
surface  from  the  drop  down  menu  and  click  on  'Add  Control  Surface'. 


nOURE  12.  CONTROL  SURFACE  GEOMETRY  SELECTION 

It  should  be  noted  from  the  above  figure,  that  there  are  a  variety  of  possible  control 
surfaces.  A  sample  sketch,  a  description  of  the  required  parameters,  and  the  generating 
equations  are  shown  below  for  each  of  the  available  control  surfaces. 
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Control  Surface  Parameters 

Before  specifying  the  geometry  of  the  control  surfaces,  a  number  of  parameters  (see 
below)  which  describe  the  location/orientation  of  the  surfaces  on  the  missile  must  be 
specified.  These  parameters,  along  with  a  description  of  each,  are  shown  below. 


HGURE  13.  CONTROL  SURFACE  PARAMETERS 

• .  Label  section  -  Enter  a  text  label  that  describes  the  control  surface. 

•  Hinge  line  -  Enter  the  axial  distance  to  the  control  surface  hinge  line  measured 
from  the  missile  nose.  The  hinge  line  is  the  point  on  the  control  surface  about 
which  deflections  are  made. 

•  Overlap  previous  control  surface  -  Select  this  option  when  two  different  types 
of  control  smfaces  are  desired  at  a  given  axial  location  (i.e.,  different  size 
vertical  and  horizontal  stabilizers).  A  detailed  description  of  this  feature  is 
provided  in  the  Linked  Fins  section. 

•  Offset  angles  -  The  offset  angles  describe  the  orientation  of  the  fins.  For  the  '  + ' 
configuration  the  corresponding  angles  are  90°,  180  ,  270,  and  360. 
Similarly,  for  the  'x'  configuration  the  angles  are  45°,  135°,  225°,  and  315 . 

Note:  1)  The  Zeus"^  code  does  not  require  any  form  of  symmetry,  and 
therefore,  any  desired  combination  of  offset  and/or  deflection  angles 
may  be  specified. 
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2)  The  control  surfaces  must  lie  on  a  zonal  boundary.  Therefore, 
the  zone  orientations  are  equivalent  to  die  offset  angles.  In  other 
words,  the  zone  edges  will  occur  at  the  above-specified  offset 
angles. 

•  Deflection  angles  -  The  deflection  angles  are  specified  in  a  clockwise  manner. 
In  other  words,  specifying  +10  for  both  the  90°  and  270°  fins  would  pitch  the 
former  up  and  the  latter  down  10°.  If  a  symmetric  deflection  is  desired,  the 
required  inputs  are  +10  for  the  90°  and  -10  for  the  270°  fin. 

Note:  The  Zeus^  code  does  not  require  any  form  of  symmetry,  and  therefore, 
any  desired  combination  of  offset  and/or  deflection  angles  may  be 
specified. 

•  Aileron/Slat  Angles  -  In  addition  to  deflecting  the  entire  surface  about  the  hinge 
line  (see  above),  it  is  also  possible  to  deflect  a  portion  of  a  control  surface.  This 
feature  allows  for  the  modeling  of  either  Aileron  or  Slat  type  controls.  To 
model  either  a  slat  or  an  aileron,  first  select  on  the  'Model  Slats/ Ailerons'  and 
then  enter  the  appropriate  deflections  in  the  'Ailerons/Slats  Angles'  dialog  box. 
Finally,  enter  &e  'ilocation'  of  either  the  slat  termination  point  or  the  aileron 
initiation  point.  The  'iloc'  parameter  corresponds  to  the  index  of  the 
terminatioiifinitiation  point  along  each  plane  defining  the  control  surface  and  is 
explained  in  the  User  Defined  Geometry  section. 

Note:  When  modeling  either  flaps  or  slats,  it  is  not  possible  to  deflect  a  portion 
of  the  trailing/leading  edges.  The  entire  edge  is  deflected  about  the  'iloc' 
point  by  the  amount  specified  in  the  Aileron/Slat  Angles  dialog  box. 

•  Delete  a  fin  -  The  number  of  fins  is  equal  to  the  maximum  number  of  zones 
specified  in  the  grid  parameter  section.  However,  it  is  possible  to  delete  a  fin 
by  right-mouse  clicking  on  the  desired  fin.  Deletion  causes  the  offset, 

deflection,  and  the  Aileron/Slat  angles  to  be  replaced  with ' - ’  to  indicate  that 

they  have  been  eliminated.  To  add  the  fin  back  into  the  calculation,  simply  right 
click  on  it  a  second  time.  This  feature  is  useful  if  a  variable  number  of  fins  are 
desired  for  the  different  sets  of  control  surfaces  (i.e.,  2  canards  and  4  tail  fins). 


3.3.3  User  Defined  Geomet 


The  most  general  method  of  specifying  the  control  surface  geometry  is  the  'User 
Defined  Geometry'  option.  For  this  method,  the  actual  geometry  points  (x,  y,  z)  are 
prescribed  by  the  user.  The  generality  of  this  method  allows  nearly  any  form  of  synunetric 
control  surface  to  be  examined.  Two  example  cases  will  be  discussed  in  the  current 
section,  the  first  of  which  is  a  simple  modified  wedge  control  surface.  The  figure  shown 
below  represents  the  upper  half  of  the  symmetric  control  surface  and  would  mate  to  the 
missile  body  along  edge  '1-2-3-4'. 
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nGURE14.  SAMPLE  1— USER  DEFINED  CONTROL  SURFACE  PARAMETERS 

The  Zeus"^  tool  requires  that  the  control  smfaces  be  described  by  a  series  of  planar 
surfaces  (see  'A',  ’B’,  and  'C  above).  In  other  words,  'A',  'B',  and  'C  each  have  a 
single  unit  outer  normal.  This  allows  the  code  to  calculate  unit  normals  to  each  section,  and 
hence,  to  generate  a  computational  domain  around  the  control  surface.  The  geometry  is 
specified  by  decomposition  into  a  series  of  planes,  all  of  which  must  have  the  same  number 
of  node  points.  For  the  figure  shown  above  the  surfaces  'A',  'B',  and  'C  are  described 
by  specifying  two  planes  with  four  points  on  each  plane  (Points  '1',  '2',  '3',  and  '4' 
define  Plane  #1,  Points  '5',  '6',  7’,  and  '8'  define  Plane  #2).  The  first  plane  (Plane  #1)  is 
that  which  corresponds  to  the  root  chord  (i.e.,  lies  along  the  missile  surface).  Successive 
planes  march  outwards  until  the  final  plane  (Plane  #2)  is  reach  at. the  tip  chord. 

Note:  1)  The  coordinate  system  for  entering  the  control  surface  points  is  based  on 
a  fiin  local  system. 

'x'  is  the  thickness  of  the  control  surface. 

'y'  is  the  radial  distance  measured  relative  to  the  missile  surface. 
(y=0  corresponds  to  the  surface  of  the  missile). 

'z'  is  the  axial  distance  measured  relative  to  the  hinge  line. 

2)  The  control  surfaces  must  be  symmetric,  and  therefore,  only  the  upper 
half  of  the  control  surface  is  specified  in  the  'User  Defined  Geometry' 
section.  Because  of  this  symmetry  requirement,  the  leading  and  trailing 
edges  of  the  fin  must  have  zero  thickness  (x=0).  This  requires  that  x  be 
zero  for  the  first  and  last  points  on  each  plane  describing  the  control  surface. 

3)  Each  plane  must  contain  the  same  number  of  node  points  (four  for  the 
above  example). 

The  above  example  can  be  used  to  model  either  slats  or  flaps.  Note  that  four  points 
are  used  to  define  each  of  the  planes  (Plane  #1  and  Plane  #2).  To  model  slats,  set  the  'iloc' 
parameter  to  '2',  select  the  'Fore  of  iloc  (Slats)'  dialog,  and  specify  the  slat  deflection  angle 
in  the  Control  Surface  Parameters  section.  This  causes  all  points  ahead  of  'iloc'  (surface 
'A')  to  be  deflected  by  the  specified  slat  angle.  Similarly,  to  model  flaps,  set  'iloc'  equal  to 
3  and  select  the  'Aft  of  iloc  (Flaps)'  dialog  to  deflect  all  points  aft  of  the  third  planar  point 
(i.e.,  deflect  surface  'C'). 

The  second  example  represents  a  more  complex  control  surface  where  three  planes 
must  be  considered  to  properly  define  the  geometry. 


19 


NSWCDD/TR-98/147 


Plane  #l 


Plane  #2 


Plane  #3 


HGURE  15.  SAMPLE  2— USER  DEFINED  CONTROL  SURFACE  PARAMETERS 


From  the  above  figure  it  is  seen  that  the  surfaces  ’  A/  and  'A2'  are  not  equivalent 
planar  surfaces  (i.e.,  do  not  have  the  same  unit  outer  normal).  A  similar  observation  is 
made  ofthe  trailing  edge  of  the  fin  (surfaces  'Bj'  and  'Bj').  Therefore,  the  above  surface 
cannot  be  described  using  only  Plane  #1  and  Plane  #3  shown  above.  Plane  #2  must  be 
included  to  properly  describe  the  control  surface. 

Note:  Each  plane  must  contain  the  same  number  of  nodal  points.  Plane  #1  contains 
four  nodal  points;  therefore,  Planes  #2  and  #3  must  also  contain  four  points.  To 
meet  this  constraint,  simply  enter  the  center  point  for  planes  '  2 '  and  ’  3 '  twice  (i.e., 
duplicate  points  '5'  and  '7')- 

Entry  of  the  'User  Defined  Geometry'  is  performed  in  the  following  dialog  box. 
First,  enter  all  of  the  'x'  points  for  a  given  plane,  then  aU  of  the  'y ' ,  and  finally,  all  of  the 
'z'  points.  Click  the  'Next  Plane'  button  to  continue  the  process  until  all  planes  have  been 
entered.  Click  on  'Finished'  after  entering  all  of  the  geometry  points. 


HGURE  16.  USER  DEFINED  CONTROL  SURFACE  INPUTS 


•  'x 'y 'z '  points  -  Local  fin  coordinates  used  for  specifying  the  control 
surface  geometry,  'x'  is  the  fin  thickness,  'y'  is  the  radial  distance  measured 
fi-om  the  missile  surface,  and  'z'  is  the  axial  distance  measured  relative  to  the 
hinge  line. 
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•  Planes  -  Displays  the  current  plane  as  well  as  the  maximum  plane  number. 


3.3.4  Modified  Wedge 


If  the  control  surface  is  a  modified  wedge,  this  section  can  be  used  in  place  of 
the’User  Defined  Geometry’  to  simplify  the  geometric  inputs  required  by  the  user. 


FIGURE  17.  MODIFIED  WEDGE  CONTROL  SURFACE  PARAMETERS 

•  Half  Angles  -  Enter  both  the  leading  and  trailing  edge  half  angles  (Uj  and  aj, 
respectively). 

•  Sweep  Angle  -  Enter  the  leading  edge  sweep  angle  (b) 

•  Span  -  Enter  the  control  surface  span  (b) 

•  Distance  to  hinge  line  -  Enter  the  distance  (on  the  root  chord)  from  the  leading 
edge  to  the  hinge  line. 

•  Wedge  thickness  -  Enter  the  leading/trailing  edge  wedge  thickness  at  both  the 
root  and  the  tip  of  the  control  surface  (tj  and  tj,  respectively). 

•  Chord  -  Enter  the  root  and  tip  chords  (Cj  and  Cj,  respectively) 

•  Convert  -  When  entering  a  control  surface,  it  is  often  desirable  to  make 
modifications  to  the  types  of  control  surfaces  allowed  in  Zeus’’^.  If  the  control 
surface  under  consideration  is  similar  to  a  modified  wedge  then  start  by  entering 
all  of  the  aforementioned  parameters,  then  press  the  'Convert  to  User  Specified 
Geometry’ button  to  convert  the  wedge  parameters  to  ’x’,  ’y’,  ’z’  coordinates. 
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It  is  then  possible  to  modify  the  individual  points  to  obtain  the  exact  control 
surface  desired. 

3.3.5  NACAOOxx 


This  section  is  used  to  describe  a  symmetric  NACAOOxx  airfoil.  Note  that  due  to 
the  bluntness  of  the  leading  edge,  difficulty  will  most  likely  be  encountered  if  this  fin  is 
evaluated  in  the  Euler  solver.  It  can,  however,  be  used  to  generate  a  Navier-Stokes  grid. 


HGUEtE  18.  NACA  4-DIGIT  AIRFOIL  PARAMETERS 


•  Number  of  points  -  Enter  the  number  of  points  used  to  describe  the  airfoil 
geometry.  Note  that  increasing  this  number  will  not  significantly  increase  the 
computational  time.  It  is  simply  the  number  of  points  used  to  draw  the  two- 
dimensional  airfoil  curve. 

•  Thickness  ratio  -  Enter  the  thickness  ratio  as  a  percentage  of  the  chord,  (i.e., 
12%  refers  to  a  thickness  of  0.12*chord) 

•  Sweep  Angle  -  Enter  the  leading  edge  sweep  angle  in  degrees. 

•  Span  -  Enter  the  span  (i.e.,  distance  firom  root  to  tip). 

•  Root/tip  chord  -  Enter  both  the  root  and  tip  chords. 

•  Distance  to  hinge  line  -  Enter  the  normalized  distance  (h/c)  fi:om  the  leading 
edge  to  the  hinge  line  (0-1)  along  the  root  chord. 

•  Convert  -  When  entering  a  control  surface,  it  is  often  desirable  to  make 
modifications  to  the  types  of  control  surfaces  allowed  in  Zeus"^.  If  the  control 
surface  under  consideration  is  similar  to  a  NACA  4-digit  airfoil  then  start  by 
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entering  all  of  the  aforementioned  parameters,  then  press  the  'Convert  to  User 
Specified  Geometry’  button  to  convert  the  fin  parameters  to  a  set  of  ’  x ' ,  ’  y ' , '  z  ’ 
coordinates.  It  is  then  possible  to  modify  the  node  points  to  obtain  the  exact 
control  surface  desired. 


The  following  is  the  equation  used  to  calculate  the  symmetric  NACAOOxx  airfoils 
(where  xx  refers  to  the  maximum  thickness  ratio): 

f  =  t(a„z®-5-a,z-a,z2+a,z3-a,z^) 

ao=  1.4845,  ai=  0.6300 
a2  =  1.7580,  a3=  1.4215 
a4  =  0.5075 
_  _  r  _  _  z 
c’  c 

'r'  is  the  local  radius,  'z'  is  the  axial  distance,  'c'  is  the  chord,  and  'f  is  the  airfoil 
thickness  ratio  (thickness/chord). 

3-3.6  Linked  Fins 

It  is  often  desirable  to  consider  different  types  of  control  surfaces  at  a  given  axid 
location  (i.e.,  different  size  vertical  and  horizontal  stabilizers).  Zeus"^  implements  ^s 
feature  by  linVing  two  sets  of  control  surfaces  together,  as  described  by  the  following 
example. 

1.  Define  the  first  control  surface.  For  example  purposes  create  a  4  finned 
configuration  oriented  in  the  '  +  '  configuration  (i.e.,  fins  at  90°,  180°,  270°,  and 
360°). 

2 .  Delete  the  90°  and  270°  fins  on  the  first  control  surface. 

3 .  Define  a  second  control  surface. 

4.  Delete  the  opposing  fins  on  the  second  control  surface  (i.e.,  the  180°  and  360°). 

5.  Select  the  “Overlap  Previous  Control  Surface”  dialog  button  for  the  second  set  of 
control  surfaces. 
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3.4  GRID  GENERATION 
Menu  Options 

In  this  section,  two  separate  computational  domains  are  examined:  the  Euler  grid 
and  the  Navier-Stokes  grid.  The  Euler  grid  is  the  one  used  by  the  Zeus‘S  tool  (see  notes 
below)  for  the  inviscid  and  boundary  layer  solvers.  The  Navier-Stokes  grid  is  generated, 
along  with  the  appropriate  boundary  condition  files;  however,  this  grid  is  not  used  by  the 
Zeus"^  tool  in  the  flowfield  analysis.  Rather,  this  grid  (and  be  file)  is  created  for  use  in  a 
full  Navier-Stokes  solver.  A  number  of  menu  options  are  available  and  are  explained 
below.  The  first  deals  with  the  inviscid  Euler  grid  and  the  remaining  three  pertain  to  the 
Navier-Stokes  grid. 

Euler  Grid 
Generate  N-S  Grid 


Normal  Spacing 
View  N-S  Grid 
3.4.1  Euler  Grid 

3. 4. 1.1  View  Grid  ('Optional!.  Before  running  the  inviscid  code,  it  is  often 
desirable  to  examine  the  computational  domain  to  ensure  that  the  relevant  geometric  features 
are  spatially  resolved  and  to  ensure  that  the  control  surfaces  are  input  correctly.  Two 
methods  of  viewing  the  computational  domain  are  available  and  are  shown  in  the  dialog 
box  below.  In  the  first,  the  grid  is  viewed  as  a  series  of  two-dimensional  planes  (see 
sample  grid  below).  The  second  method  displays  the  Euler  grid  in  three-dimensional  block 
format. 


Generate  and  view  the  Euler  grid 

Generate  the  Navier-Stokes  grid  and  boundary  condition 

data  files  'filename\RunGrid\GASPGrid.p3da'  & 

'GASPBC.inp' 

View  die  radial  spacing  (at  the  first  point  off  the  wall) 

along  the  missile 

Plot  the  Navier-Stokes  grid. 


When  using  the  '2d  planes'  option,  it  is  possible  to  generate  a  movie  of  the  grid  that 
steps  axially  along  the  missile.  To  generate  the  movie  file  from  the  two-dimensional  grid, 
run  the  Tecplot  macro  'Playback\2dGrid_Movie.mcr'.  The  movie  file  is  written  to  the  file 
'Playback\Zeus_Grid2d.rm'  and  can  be  replayed  using  the  Tecplot  Framer  program. 

Note:  1)  The  axial  stepsize  taken  during  the  inviscid  computation  is 
controlled  by  the  Courant-Friederichs-Lewy  (CFL)  stability  condition 
specified  in  the  integration  parameter  section.  This  is  calculated  at  run-time, 
and  therefore,  the  axial  steps  taken  during  the  grid  generation  are  not 
equivalent  to  those  taken  during  the  actual  solution.  For  the  purpose  of 
evaluating  the  computational  domain,  equal  axial  spacing  is  taken  along  the 
missile. 

2)  There  are  two  possible  methods  for  generating  the  outer  boundary 
of  the  computational  domain.  In  the  first,  it  is  specified  in  terms  of  body 
radii  in  the  grid  parameter  section.  In  the  second,  the  outer  boundary  tracks 
the  conical  shock  as  the  solution  proceeds  downstream.  For  the  latter 
method,  the  conical  shock  location,  and  thus  the  outer  boundary  is  not 
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known  until  rantime,  and  therefore,  for  the  purposes  of  grid  generation,  the 
first  method  of  calculating  the  outer  boundary  will  be  employed. 

3)  Viewing  the  Euler  grid  is  not  required  before  executing  the 
inviscid  solution.  The  actual  grid  used  during  the  inviscid  computation  is 
calrnl^tpd  at  runtime  and  can  be  viewed  using  the  post-processing  features 
of  the  Zeus"^  code. 


nGURE19.  EULER  GRID  GENERATION  PARAMETERS 

•  Plot  type  -  Generate  the  grid  as  either  a  series  of  two-dimensional  planes  (i.e., 
each  plane  is  written  as  a  separate  Tecplo^  zone),  or  a  one  three-dimensional 
grid  (i.e.,  block  format). 

•  Starting/Jinal  zeta  -  Enter  the  starting  and  final  axial  locations  for  the  grid 
generation.  Any  subsection  can  be  viewed  using  the  grid  generator,  and 
therefore,  the  z^  values  do  not  have  to  correspond  to  the  starting/final  axial 
locations  of  the  missile. 

•  Number  of  steps  -  Enter  the  number  of  axial  steps  to  be  taken.  The  axial 
stepsize  is  assumed  constant  for  the  purposes  of  grid  generation  and  is 
calculated  linearly  using  the  starting  and  final  zeta  values  specified  above. 

3.4.2  Navier-Stokes  Grid 

3.4.2. 1  Generate  Grid.  The  Zeus"^  tool  can  be  used  to  generate  solutions  for 
missile  type  geometries  using  both  an  Euler  and  a  boundary  layer  solver.  However,  it  is 
often  desirable  to  perform  a  more  in-depth  analysis  using  a  full  Navier-Stokes  solver.  For 
these  cases,  the  Zeus"*^  tool  serves  as  a  grid  generator,  allowing  the  user  to  create  complex 
three-dimensional  grids  in  a  matter  of  minutes.  Currently,  the  tool  is  setup  to  generate  CFD 
grids  in  Plot3d  format  ('filename\RunGrid\GASPGrid.p3da’)  along  with  the  appropriate 
zonal  boundary  conditions  ('filename\RunGrid\GASPBC.inp')  for  the  GASP'  flow  solver. 
The  boundary  condition  file  specifies  each  point  on  a  zone  edge  as  either  a  solid  surface  (a 
fin)  or  a  zon^  boundary. 
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HGURE  20.  NAVIER-STOKES  GRID  GENERATION  PARAMETERS 

•  Axial  clustering  -  In  this  section  the  axial  clustering  for  the  Navier-Stokes  grid 
is  specified.  The  radial  and  circunrferential  values  are  declared  in  the  grid 
parameter  section. 

Note:  1)  Zeus‘S  uses  a  Newton  method  to  calculate  the  axial  clustering 
parameters  (i.e.,  geometric  stretching  factors).  To  prevent  a 
divergent  solution,  a  valid  initial  spacing  must  be  provided.  This  is 
accomplished  by  first  specifying  the  number  of  axial  sections  along 
the  missile.  Next,  enter  the  number  of  points  and  the  ending  axial 
location  for  each  section.  Fmally,  click  on  the  'Reset  Spacing' 
button  to  calculate  the  appropriate  parameters  for  an  equally  spaced 
grid.  Once  a  valid  grid  has  been  obtained,  the  beginning/ending 
spacing  can  be  adjusted  incrementally  to  obtain  the  desired  axial 
clustering. 

2)  The  grid  and  boimdary  condition  output  files  produced  during  the 
grid  generation  process  are  often  times  very  large,  and  therefore, 
may  take  some  time  to  generate.  Also,  the  process  of  determining 
the  appropriate  axial  clustering  is  an  iterative  procedure,  and  is 
independent  of  the  circumferential  layout  of  the  grid.  Therefore,  to 
minimize  the  amount  of  calculation  required  per  iteration,  set  the 
number  of  circumferential  points  to  unity  in  the  Navier-Stokes  grid 
section.  After  the  desired  axial  clustering  is  obtained,  reset  the 
number  of  circumferential  points  to  the  desired  value,  and  regenerate 
the  final  Navier-Stokes  grid. 


26 


NSWCDD/TR-98/147 


•  Number  of  axial  stations  -  Specify  the  number  of  axial  sections  in  which 
clustering  is  desired.  For  the  sample  grid  shown  below  the  six  sections 
designated  are  distinguished  by  the  vertical  tic  marks.  The  division  shown 
in  the  example  allows  for  axial  clustering  at  the  beginning/end  of  the  nose 
section,  leading/trailing  edges  of  both  sets  of  fins,  as  well  as  clustering  at 
the  aft  end  of  the  missile. 

•  Number  of  'z'  stations  -  The  number  of  axial  points  per  section. 

•  Ending  'z'  location  -  Ending  axial  location  for  each  section  [ft/m]. 

•  Reset  spacing  -  Reset  the  beginning/ending  axial  spacing  in  each  section  to 
obtain  equal  spacing  in  all  sections. 

•  Beginning  spacing  -  Beginning  axial  spacing  (Az)  for  each  section  [ft/mj. 

•  Ending  spacing  -  Ending  axial  spacing  (Az)  for  each  section  [ft/m]. 

•  Aft  grid  -  If  a  grid  aft  of  the  missile  is  desired  (subsonic  flow)  then  this 
checkbox  must  be  selected. 

•  Number  of  radial  points  -  This  refers  to  the  points  which  lie  rachally 
between  the  missile  centerline  and  the  missile  surface  directly  behind  the  missile 
(see  yellow  grid  section  or  area  designated  by  the  horizontal  tic  marks).  The 
radial  spacing  at  the  first  point  is  taken  to  be  the  same  as  that  on  the  surface  (at 
the  aft  end)  of  the  missile.  Zeus"^  calculates  the  geometric  stretching  parameter 
that  places  the  final  radial  point  along  the  missile  centerline.  If  the  error 
message  'Unattainable  radial  clustering'  is  displayed,  try  reducing  the  number 
of  radial  points  in  the  aft  grid  section. 

•  Outer  boundary  clustering  -  Specify  the  outer  boundary  clustering  around 
the  nose  section  of  the  missile. 

•  First  point  angle  -  The  angle  between  the  missile  centerline  and  the  first  nose 
point  [degrees].  See  angle  'a'  below.  The  angle  must  be  greater  than  0 
degrees. 


•  Last  point  angle  -  The  angle  between  the  missile  centerline  and  the  final  nose 
point  [degrees].  See  angle  'b'  below.  The  angle  must  be  less  than  90 
degrees. 


FIGURE  21.  SAMPLE  CROSS-SECTION  OF  NAVIER-STOKES  GRID 
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The  final  step  in  generating  a  Navier-Stokes  grid  is  to  specify  the  surface  boundary 
conditions  (see  figure  below). 


HGURE  22.  NAVIER-STOKES  BOUNDARY  CONDITION  PARAMETERS 


•  Surface  boundary  condition  -  Specify  the  type  of  boundary  condition  to  be 
employed  when  a  solid  surface  is  encountered.  Currently,  the  Zeus"^  tool 
generates  a  boundary  condition  file  for  the  GASP  flow  solver,  which 
specifies  the  boundary  condition  along  the  zone  edges  (i.e.,  solid  surface  or 
zonal  boundary). 

•  Two-Equation  Model  -  The  type  of  two-equation  turbulence  model  applied 
during  the  Navier-Stokes  solution  must  be  specified  in  the  GASP  boundary 
condition  input  file. 

3.4.2.2  View  Normal  Spacing.  A  prerequisite  for  accurately  computing  a  Navier- 
Stokes  solution  is  the  spatial  resolution  of  Ae  boundary  layer,  and  therefore,  the 
computational  grid  must  be  appropriately  refined  near  the  surface.  The  boundary  layer 
thickness  cannot  be  determined  a  priori,  however,  an  approximate  value  can  be  obtained 
using  flat  plate  assumptions.  Using  this  estimate  for  the  boundary  layer  thickness,  the 
radial  clustering  parameter  can  be  increased/decreased  until  an  appropriate  value  for  the 
normal  spacing  to  the  first  point  off  the  surface  is  obtained.  A  typical  value  for  the  required 
normal  spacing  is  l/20th  to  l/30th  of  the  approximate  boundary  layer  thickness  (i.e., 
approximately  20-30  points  in  the  boundary  layer). 

3. 4. 3. 3  View  Grid.  This  option  plots  the  Navier-Stokes  grid  in  three-dimensional 
block  format  using  the  Tecplol®  plotting  package.  A  two-dimensional  slice  of  the  three- 
dimensional  grid  can  be  obtained  by  executing  the  Tecplot  macro 
'Playback\GaspGrid.mcr'.  This  is  the  view  displayed  in  the  Generate  Gasp  Grid  section. 
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and  is  useful  for  evaluating  the  axial  clustering  parameters  as  well  as  the  spatial  resolution 
of  the  control  surfaces. 

Note:  The  control  surfaces  are  represented  by  the  contour  plot  variable 
'onfin'.  Unity  and  zero  represent  a  solid  surface  and  a  zonal  boundary, 
respectively.  To  evaluate  the  spatial  resolution  of  the  control  surfaces, 
examine  a  contour  plot  of  the  'onfm'  parameter. 
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CHAPTER  4:  INPUT  PARAMETERS 


OVERVIEW 

A  variety  of  parameters  are  required  before  execution  of  the  Euler  or  boundary  layer 
solvers  can  proceed.  The  runtime  parameters  are  specified  in  the  'Inputs'  section  from  the 
main  menu.  The  parameters  are  separated  into  four  basic  categories  as  shown  below. 

1 .  Aerodynamic  Data 

Run  Matrix 
Genfiral  Data 

2.  Integration  Controls 

3.  Separation  Parameters 

4.  Boundary  Layer  Controls 

Execution  Parameters 
Modeling  Parameters 
Output  Parameters 

5.  Output  Controls 


4.1  AERODYNAMIC  PARAMETERS 

4.1.1  Run  Matrix  Parameters 

The  Aerodynamic  run  matrix  is  where  the  Mach  numbers  and  missile  orientations 
are  specified.  This  is  accomplished  by  entering  an  initial,  final,  and  incremental  value  for 
Mach  number,  angle  of  attack,  angle  of  yaw,  and  roll  angle  (see  figure  below).  Zeus^  then 
calculates  an  Aerodynamic  run  matrix  consisting  of  all  possible  cases. 
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HGURE  23.  AERODYNAMIC  SWEEP  PARAMETERS 

•  Mach  number  sweep  -  Enter  the  starting,  final,  and  incremental  Mach  number. 
For  the  sample  shown  above,  a  Mach  number  of  4.0,  4.25,  and  4.50  will  be 
considered. 

•  Angle  of  attack  sweep  -  Enter  the  starting,  final,  and  incremental  angle  of  attack. 

•  Angle  of  yaw  sweep  -  Enter  the  starting,  final,  and  incremental  yaw  angle. 

•  Roll  angle  sweep  -  Enter  the  starting,  final,  and  incremental  roll  angle. 

•  Deflection  sweep  -  Running  sweeps  on  control  surface  deflections  has  not  been 
implemented  in  the  current  version. 

Note:  For  each  case  (i.e.,  every  Mach  number,  angle  of  attack,  angle  of 

yaw,  and  roll  angle  combination)  a  subdirectory  wiU  be  created  under  the  main 
directory  and  will  be  labeled  RunOOOl-Runxxxx.  In  addition,  a  file  called 
RunMatrix.txt  is  generated  which  maps  the  case  number  to  the  aerodynamic 
parameters  entered  above. 
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4. 1 .2  General  Aerodynamic  Parameters 

A  number  of  general  aerodynamic  parameters  must  be  specified  before  executing 
either  the  Euler  or  boundary  layer  solvers.  These  parameters  are  shown  in  the  following 
figure  and  described  in  detail  below. 


HGURE  24.  GENERAL  AERODYNAMIC  INPUT  PARAMETERS 

•  Free-stream  data  -  The  required  free-stream  data  includes  the  pressure,  density, 
and  ratio  of  specific  heats  (g).  There  are  three  methods  of  specifying  these 
parameters. 

1)  Select  either  of  the  standard  atmospheric  tables  (1962  U.S.  Standard 
Tables  or  1959  AFDC  Tables)  and  then  enter  an  altitude  in  the  provided 
dialog  box.  The  pressure  and  density  will  be  calculated  from  the 
standard  table  data  and  the  ratio  of  specific  heats  is  assumed  that  of 
standard  air  (1.4). 

2)  Specify  the  pressure,  density,  and  ratio  of  specific  heats  directly  by 
clicking  on  the  radio  button  next  to  these  parameters  and  entering  the 
desired  values. 
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3)  Specify  a  unit  Reynolds  number  and  let  Zeus'^  calculate  the 
corresponding  pressure  and  density.  Gamma,  the  ratio  of  specific 
heats,  is  taken  to  be  that  of  standard  air  ( 1 .4). 


•  Units  -  Specify  either  English  or  metric  units. 

Note:  All  of  the  parameters  entered  must  be  in  base  units  (i.e.,  meters,  not 
centimeters;  feet,  not  inches,  ...etc.).  Incorrect  results  will  be  obtained 
if  you  enter  all  of  your  lengths  in  centimeters  and  then  try  to 
compensate  by  entering  the  reference  quantities  in  centimeters. 
Although  entering  non-standard  units  will  not  effect  the  inviscid 
solution,  it  will  produce  incorrect  results  for  the  viscous  solution  due  to 
the  Reynolds  number  effect. 

•  Reference  conditions  -  Enter  the  reference  diameter  (typically  the  base  diameter) 
and  the  reference  area  (typically  the  area  calculated  using  the  base  of  the 
missile).  Also,  specify  the  axial  location  about  which  the  pitching  and  yawing 
moments  are  calculated  (The  rolling  moment  is  always  taken  about  the  missile 
centerline). 

Note:  A  nose  down  pitching  moment  is  negative. 

•  Viscous  parameters  -  If  the  viscous  boundary  layer  code  is  to  be  considered 
then  the  wall  boundary  condition  for  the  heat  flux  and  the  Prandtl  number  must 
be  specified.  For  the  adiabatic  wall,  calculating  an  appropriate  wall  temperature 
using  the  recovery  factor  methodology  imposes  a  zero  heat  flux  boundary 
conation. 

•  Gas  type  -  Select  either  a  perfect  gas  or  a  real  gas  solution. 

Note:  If  'real  gas'  is  selected  then  the  approximate  Riemann  solver  must  be 
selected  in  the  integration  controls  section. 


4.2  INTEGRA’nON  CONTROLS 
4.2.1  Parameter  Description 

A  number  of  user  specified  integration  parameters  must  be  specified  in  the 
following  section.  These  include  axial  limits  for  integration,  starting  solution,  stability 
criterion,  limiters,  and  flux  calculation  parameters. 
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nGURE25.  NUMERICAL  INTEGRATION  CONTROLS 

•  Space  marching  -  Enter  the  starting  and  final  axial  locations  for  the  solution. 
Also,  enter  a  maximum  number  of  steps  to  take  before  the  code  aborts. 

Note:  The  starting  axial  location  is  dependent  on  the  type  of  starting  solution 
applied  (see  below  for  an  explanation  of  each  of  &e  different  methods  of 
generating  a  starting  solution). 

•  Stability  criterion  -  Enter  a  step  size  safety  factor  (analogous  to  a  CFL  number) 
between  zero  and  unity.  Also,  specify  a  maximum  allowable  step  size  to  ensure 
that  all  relevant  features  are  spatially  resolved  in  the  axial  direction. 

•  Starting  solution  -  The  Euler  solver  is  a  marching  code,  and  therefore,  a  starting 
solution  must  be  specified  on  the  missile  before  the  integration  can  proceed. 
Three  different  types  of  starting  solutions  are  available,  each  of  which  is 
discussed  below. 

1 )  Free-stream  -  A  fi-ee-stream  profile  is  specified  as  the  starting  solution  at  the 
axial  location  entered  in  the  'space  Marching'  section.  For  this  type  of  inflow 
profile,  the  outer  boundary  does  not  track  the  conical  shock,  but  rather,  is  user 
specified  in  terms  of  body  radii  in  the  grid  parameter  section. 

Note:  The  region  upstream  of  this  profile  is  ignored,  and  therefore,  in  order  to 
accurately  compute  the  forces  and  moments,  the  integration  should  begin  as 
close  to  Ae  nose  as  is  numerically  possible. 
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2)  Cone  Solution  -  The  cone  starting  solution  applies  a  one-dimensional  cone 
solver  to  the  region  from  the  nose  to  the  beginning  axial  location  specified  in  the 
space  marching  section.  This  calculates  the  shock  angle,  as  well  as  the 
flowfield  directly  aft  of  the  shock,  given  the  local  turning  angle.  The  starting 
profile  is  then  taken  to  be  that  generated  by  the  conical  solver.  In  other  words, 
a  conical  flow  solver  is  apphed  from  the  nose  tip  to  the  location  specified  as  the 
'Axial  Location  to  begin  computation'.  The  Euler  solver  then  begins  integration 
at  this  location  with  the  profile  supplied  by  the  conical  flow  solver. 

Note:  a)  The  conical  flow  solver  is  intended  to  provide  the  shock  angle  as  well 
as  the  flowfield  properties  directly  aft  of  the  shock.  It  is  not  intended  to 
solve  the  flowfield  for  the  entire  nose  cone.  The  one-dimensional  cone 
solver  cannot  calculate  the  forces/moments,  and  therefore,  if  you  do  not 
start  the  Euler  solver  until  the  end  of  the  nose  section,  the 
forces/moments  due  to  this  section  will  be  neglected.  In  other  words,  the 
conical  solver  should  only  be  applied  to  a  very  small  region  around  the 
nose  tip  and  not  to  the  entire  nose  section. 

b)  If  the  local  turning  angle  becomes  excessively  large  then  the  one¬ 
dimensional  cone  solver  will  compute  a  subsonic  Mach  number  aft  of  the 
shockwave.  The  Zeus  Euler  solver  is  a  space  marching  code,  and 
therefore,  cannot  handle  subsonic  regions.  For  this  case,  a  free-stream 
starting  solution  must  be  applied. 

3)  Blunt  body  -  For  a  blunt  body,  the  flow  will  be  subsonic  in  some  region 
surrounding  the  nose  section.  If  you  attempt  to  begin  the  Zeus  integration  too 
close  to  the  nose,  the  code  will  encounter  a  subsonic  region  and  abort.  To 
prevent  this,  a  blunt  body  solver  must  be  applied  to  the  entire  subsonic  region. 
This  solver  provides  a  solution  from  the  beginning  of  the  nose  to  the  value 
specified  in  the  space  marching  section  labeled  'Axial  location  to  be^ 
computation'.  The  user  must  ensure  that  this  value  is  downstream  of  the  entire 
subsonic  region. 

Note:  In  contrast  to  the  free-stream  and  conical  starting  profiles,  the  blunt 
body  solver  calculates  the  forces  and  moments  in  the  region  between  the 
nose  and  the  starting  location  for  the  Euler  solver.  These  forces  and 
mnrppntR  are  summed  with  the  values  obtained  by  the  Euler  solver  to 
provide  the  total  forces/moments  acting  on  the  body.  Thus,  for  this 
starting  procedure,  it  is  acceptable  to  run  the  blunt  body  solver  over  the 
entire  nose  section. 

•  Riemann  solver  -  Specify  whether  the  fiill  Riemann  or  approximate  Riemann 
solver  is  applied. 

•  Bow  shock  flux  -  Calculate  the  numerical  flux  at  the  bow  shock  using  either  the 
complete  Riemann  solver  or  freestream  properties. 
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•  Bow  shock  angle  -  Calculate  the  bow  shock  angle  using  either  the  complete 
Riemann  solver  or  the  Roe  averaged  variables. 

•  Interior  point  limiter  -  Numerical  limiter  at  all  interior  points  (i.e.,  not  adjacent 
to  a  surface).  The  default  value  is  unity,  however,  this  can  be 
increased/decreased  to  modify  the  numerical  damping  provided  by  the 
integration  scheme.  The  applicable  range  is  2  to  0,  where  2  and  0  are  the  least 
and  most  dissipative,  respectively.  A  value  of  zero  reduces  the  integration 
scheme  to  a  first-order  Gudonov  method. 

Note:  Decreasing  the  damping  (i.e.,  increasing  the  limiter  value  above)  will 
improve  the  accuracy  of  the  numerical  scheme  and  will  more  accurately 
cq)ture  the  flowfield  discontinuities.  Insufficient  damping  can  often 
lead  to  numerical  oscillations  in  the  flowfield  and  prevent  a  stable 
integration. 

•  Turning  angle  correction  -  If  a  subsonic  region  is  encountered  while  turning  the 
flow,  the  Zeus  code  will  apply  a  limiter  known  as  a  turning  angle  correction. 
The  actual  turning  angle  encountered  in  the  flow  is  multiplied  by  the  turning 
angle  correction  to  try  to  prevent  a  subsonic  region.  If  this  correction  does  not 
prevent  a  subsonic  region  then  the  code  will  abort.  Typically,  a  value  of  0.9 
(i.e.,  90%  of  the  original  angle)  is  applied. 


4.3  SEPARATION  PARAMETERS 
4.3.1  Parameter  Description 

The  Zeus  code  is  an  inviscid  solver,  and  as  such,  neglects  the  effects  of  viscosity 
on  die  flowfield.  Since  flow  separation  is  a  viscous  phenomenon,  an  Euler  solver  alone 
cannot  accurately  predict  the  effects  of  a  separation  zone  on  the  flowfield.  Two  separation 
models  (clipping  and  forced)  are  available  in  the  Zeus  Euler  solver  which  attempt  to 
improve  the  predictive  capabilities  of  the  tool  near  separation  zones. 

The  clipping  separation  model  is  described  in  detail  in  Ref.  1 1  and  operates  by 
decreasing  the  crossflow  velocity  on  and  near  the  body  surface.  The  crossflow  velocity 
reduction  is  accomplished  by  setting  an  upper  limit  to  the  allowable  crossflow  velocity.  If 
the  velocity  at  any  point  exceeds  this  value,  it  is  reduced  to  this  level.  Pressure  and  density 
are  assumed  imchanged  and  the  axial  component  of  velocity  is  increased  to  give  the  correct 
stagnation  enthalpy  value.  Clipping  destroys  the  crossflow  shock  and  produces  a  large 
vortex  on  the  leeside  of  the  body  that  is  in  qualitative  agreement  with  experiment. 
Computed  pressures  on  the  leeside  of  the  body  are  in  better  agreement  with  experiment. 
However,  clipping  tends  to  increase  the  windward  pressures  and  often  decreases  the 
accuracy  of  the  results  in  this  region. 

The  forced  separation  model,  originally  described  in  this  form  by  Ref.  12,  seeks  to 
simulate  separation  by  altering  the  velocity  direction  along  a  user  defined  separation  line. 
In  each  crossflow  plane,  the  separation  model  is  applied  to  wall  cells  “s”  and  “s+l”,  which 
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are  inunediately  windward  and  leeward  of  the  separation  line,  respectively.  Pressure  and 
density  at  these  points  are  defined  by: 


Ps  =^-^(p'sA+ P's^x) 
Ps  =  0.5(p'_i  +  Pj+i) 

PsA  =^-^(p's  +  Pm) 


Here '  denotes  old  values.  The  flow  velocities  are  determined  by  prescribing  a  streainline 
direction  and  assuming  constant  stagnation  enthalpy.  Unfortunately,  solutions  obtained 
with  the  forced  separation  model  are  sensitive  to  mesh  size.  As  the  mesh  is  refined,  the 
pressure  beneath  the  vortex  diminishes. 


FIGURE  26.  SEPARATION  MODELING  PARAMETERS 

•  Separation  Model  -  Select  the  desired  separation  model. 

•  Clipping  Separation  -  Enter  the  axial  location  to  begin  clipping. 

•  Forced  Separation  -  Enter  the  number  of  forced  separation  lines.  For  each 
of  the  lines,  the  user  must  specify  the  following  parameters: 

•  ISIDE  -  The  side  of  the  body  on  which  the  separation  occurs.  Zero  if 
separation  occurs  between  0°  and  180°,  unity  if  separation  occurs  between  180° 
and  360°. 
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•  ISEP  -  The  number  of  points  used  to  define  the  separation  line. 

•  ZSSEP  -  Zeta  value  at  which  separation  is  started. 

•  ZESEP  -  Zeta  value  at  which  separation  is  terminated. 

•  PHICD,  PHIAD,  BETACD,  BET  AD  -  Flow  direction  in  degrees.  Typical 
values  are  20°,  20°,  20°,  5°,  respectively. 

•  ZEPZyZEPP  -  A  list  of  ISEP  pairs  of  separation  line  coordinates  where 
ZEPZ  and  ZEPP  are  the  z  and  (p  coordinates  for  each  point. 


4.4  BOUNDARY  LAYER  CONTROLS 
4.4. 1  BL  Execution  Parameters 

The  boundary  layer  code  requires,  as  a  boundary  condition,  the  inviscid  solution 
from  the  Zeus  code.  Therefore,  before  executing  the  boundary  layer  code,  Zeus^  must 
save  the  inviscid  solution  at  a  number  of  different  axial  stations.  It  is  both  unnecessary  and 
infeasible  (due  to  computational  constraints)  to  save  the  inviscid  flowfield  at  every  axial 
station,  however,  an  adequate  number  of  axial  stations  must  be  considered  in  order  to 
accurately  solve  the  boundary  layer  equations.  The  user  must  save  the  inviscid  solution 
often  enough  to  resolve  all  relevant  geometric  features  axially  (i.e.,  must  have  several  axial 
stations  written  out  on  each  control  surface). 

•  OMipwrinvisc/dyZovv/ieW- Enter  the  increment  and  the  axial  stepsize  for  saving 
the  inviscid  flowfield.  The  Zeus  code  will  save  the  inviscid  flowfield  every 
incremental  axial  step  as  well  as  every  time  the  axial  stepsize  exceeds  the 
specified  value. 

Note:  The  Zeus"^  code  must  save  the  inviscid  flowfield  often  enough  to 
axiaUy  resolve  aU  relevant  geometric  features. 

•  Starting  axial  location  for  integration  -  Enter  the  axial  location  to  begin  the 
boundary  layer  solver.  Note  that  this  must  lie  downstream  of  the  starting 
location  of  the  Euler  solver  as  specified  in  the  Integration  Parameters  section. 

•  Minimum  number  of  steps  -  Enter  the  minimum  number  of  axial  steps  the 
boundary  layer  solver  will  take.  The  code  uses  this  parameter  to  calculate  a 
maximum  allowable  stepsize. 

•  Maximum  number  of  steps  -  Enter  the  maximum  number  of  axial  steps  the 
boimdaiy  layer  solver  will  take  before  aborting. 
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HGURE27.  BOUNDARY  LAYER  EXECUTION  PARAMETERS 

•  Implicit  iterations  -  Niomber  of  implicit  iterations. 

•  Effective  radius  -  This  parameter  selects  tiie  manner  in  which  the  effective 
radius  is  calculated;  fee  method  by  which  streamline  spreading  is  determined. 
The  “streamline”  option  uses  fee  distance  between  successive  streamlines,  fee 
“geometry”  option  computes  the  local  surface  radius,  while  “none”  sets  fee 
radius  to  unity.  Optimal  results  are  obtained  with  fee  “streamline”  option  except 
on  circular  bodies  at  zero  incidence.  Here  fee  “geometry”  option  is 
recommended. 
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4.4.2  BL  Modeling  Parameters 


HGURE  28.  BOUNDARY  LAYER  MODELING  PARAMETERS 

•  Transition  parameters  -  Select  one  of  the  three  possible  transition  parameters.  If 
laminar  is  selected  then  the  flow  is  forced  to  remain  laminar  for  the  entire 
calculation.  If  an  'Axial  transition  location'  is  specified,  the  flow  is  laminar  up 
to  this  point  and  then  abruptly  transitions  to  turbulent  flow.  If  the  'Gradual 
transition  factor'  is  selected,  then  the  transition  model  is  taken  from  Ref.  6,  and 
the  parameter  specified  is  F,  (default  of  1.0)  from  Eq.  23. 

•  Pressure  gradient  for  separation  -  If  the  streamwise  pressure  gradient  exceeds 
this  user  specified  parameter  the  boundary  layer  and  energy  thicknesses  are 
reduced  to  simulate  separation®. 

•  Streamline  divergence  -  Controls  streamline  divergence  as  described  in 
Reference  6  (K^  from  Eq.  32).  Set  to  unity  to  follow  inviscid  streamlines.  Use 
a  value  of  2.0  for  bodies  at  incidence.  Only  applies  to  edge  #1  of  each  zone. 

•  Streamline  convergence  -  Controls  streamline  convergence  as  described  in 
Reference  6  (K^  from  Eq.  32).  Set  to  unity  to  follow  inviscid  streamlines.  Use 
a  value  of  0.2  for  bodies  at  incidence.  Only  applies  to  edge  #1  of  each  zone. 

•  Induced  pressure  limiter  -  The  boundary  layer  code  attempts  to  correct  for  the 
induced  pressure  effects  caused  by  the  formation  of  the  boundary  layer.  To 
prevent  large  pressure  corrections  around  leading/trailing  edges,  the  induced 
pressure  is  limited  using  the  'Induced  pressure  limiter'  as  discussed  in  Ref.  6 
(Kf  in  section  4.3). 
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•  Momentum/Energy  thickness  limiter  -  A  parameter  which  limits  the 
momentum/energy  thickness  in  anticipated  regions  of  separation.  The  liimter  is 
discussed  in  Ref.  6  (k  from  Eq.  33). 

4.4.3  RL  Output  Controls 

The  following  dialog  box  controls  the  output  from  the  boundary  layer  solver. 


HGURE  29.  BOUNDARY  LAYER  OUTPUT  PARAMETERS 

•  Print  Frequency  -  The  ASCII  output  table,  'ZeusBL.out'  is  generated  every 
'print  frequency'  steps.  The  'reduced/normal/diagnostic'  option  controls  the 
amount  of  data  written  to  the  output  tape.  “Reduced”  prints  only  the  most 
important  quantities  such  as  Cf  and  Q.  “Normal”  includes  additional  boundary 
layer  parameters  as  weU  as  edge  conditions. 

•  Plot  Frequency  -  The  contour  plot  tape  is  written  every  “plot  frequency”  steps. 

•  Axial  Plots  -  Specify  the  surface  number,  point  location,  dependent  and 
independent  variables,  as  well  as  x/y  scaling  factors.  Each  of  the  axial  plots 
will  be  written  to  a  file  ‘graphxxx’,  where  ‘xxx’  is  the  number  of  the  axial  plot. 

i  -  surface  number 
j  -  point  number 


Y  variable  -  Select  one  of  the  following  for  a  dependent  variable: 
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1.  pressure 

2.  density 

3.  edge  velocity 

4.  temperature 

5.  Cf  (base  on  edge  conditions) 

6.  Cf  (base  on  reference  conditions) 

7.  Ch  (base  on  edge  conditions) 

8.  Q  (base  on  reference  conditions) 

9.  Momentum  thickness 

10.  Displacement  thickness 

11.  Energy  thickness 

12.  Heating  rate  (B/(ft'^2  sec),  J/(m''2  sec)) 

13.  Viscous  loads;  Here  “i”  is  the  surface  number  and  'j '  is  the  load 
(l=x  force,  2=y  force,  3=z  force,  4=x  moment,  5=y  moment, 
6=z  moment).  To  get  total  body  viscous  loads  set  “i”  to  0. 

14.  Induced  loads.  Same  options  apply  as  for  viscous  loads. 

Y  scale  -  Dependent  variable  scale  factor. 

X  variable  -  Select  one  of  the  following  for  an  independent  variable: 

1.  Axial  location,  z 

2 .  Reynolds  number  (based  on  axial  location),  Re^ 

3 .  Streamline  length,  s 

4.  Reynolds  number  (based  on  streamline  length),  Re^ 

X  scale  -  Independent  variable  scale  factor. 

•  Contour  Plots  -  Specify  a  list  of  target  stations  (separated  by  blank  spaces) 
where  the  contour  plots  are  desired  (the  files  are  written  to  ‘graphcxxx’,  where 
‘xxx’  refers  to  the  appropriate  contour  file  number).  Then  list  the  following 
parameters  for  each  contour  plot: 


Surface  -  Enter  the  surface  number.  The  surface  numbers  are  counted  sequentially 
along  each  of  the  four  edges  for  each  zone  (i.e..  Zone  2,  Edge3 
corresponds  to  surface  number  7). 

Variable  -  Select  one  of  the  following  for  a  dependent  variable: 

1.  pressure 

2.  density 

3.  edge  velocity 

4.  temperature 
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5.  Cf  (base  on  edge  conditions) 

6.  Cf  (base  on  reference  conditions) 

7.  Ch  (base  on  edge  conditions) 

8.  Cj,  (base  on  reference  conditions) 

9.  Momentum  thickness 

1 0.  Displacement  thickness 

1 1 .  Energy  thickness 

12.  Heating  rate  (B/(ft^2  sec),  J/(m'^2  sec)) 

y-scale  -  Dependent  variable  scale  factor, 
x-scale  -  Crossflow  variable  scale  factor. 


4.5  OUTPUT  CONTROLS 
4.5.1  Parameter  Description 

This  section  controls  which  type  of  output  the  Zeus  Euler  solver  generates,  as  well 
as  the  output  frequency.  The  first  section  in  the  figure  below  controls  the  output  for  the 
surface  pressures,  forces,  moments,  and  centers  of  pressure.  The  second  section  controls 
the  output  for  generating  contour  plots  of  the  flowfield  variables.  The  final  section  controls 
both  the  ASCn  output  and  screen  residual. 


-  - . — ■ 

FIGURE  30.  OUTPUT  CONTROL  PARAMETERS 
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•  Surface  pressures,  forces,  and  moments  -  This  section  handles  the  output 
for  the  surface  pressure,  forces,  and  moments. 

•  Print  cell  coordinates  -  Write  the  cell  coordinates  (axial  and  radial  location) 
at  each  cell  where  the  data  is  written  along  with  the  surface  pressure,  forces, 
moments,  and  centers  of  pressure. 

•  Write  surface  properties  -  Specify  both  an  incremental  value  and  a 
maximum  stepsize  for  generating  the  output  data.  The  Zeus^  code  writes  the 
output  at  every  specified  incremental  value  (axially)  as  well  as  if  the  axial 
stepsize  is  larger  than  the  specified  parameter. 

•  Contour  data  -  Contour  plots  of  a  variety  of  flowfield  parameters  (pressure, 
density,  Mach  number,  etc. . .)  can  be  generated  at  various  axial  stations. 

•  'k'  Axial  stations  -  Enter  the  increment  for  generating  the  contour  plot 
output  data.  The  code  will  write  the  output  data  at  every  user  specified 
incremental  axial  station. 

•  User  defined  stations  -  Specify  the  specific  axial  locations  where  the 
contour  data  is  desired.  This  is  useful  for  generating  contour  plots  at  specific 
locations  in  the  flowfield  (i.e.,  leading/trailing  edges  of  control  surfaces). 

•  Given  stepsize  -  The  Zeus"^  code  generates  a  contour  plot  file  every  time  the 
axial  stepsize  is  greater  than  the  specified  value. 

•  ASCII  data  -  The  ASCII  data  is  written  to  the  ‘fort.9’  data  file  in  the 
'filename'  directory. 

•  Output  residual  -  Specify  the  axial  step  increment  for  writing  out  the 
residual  data.  Note,  the  code  writes  the  residual  data  to  the  screen  as  well  as  to 
the  ‘fort.9’  data  file. 

•  Crossflow  plane  -  Specify  the  axial  increrrient  for  saving  the  ASCII  data 
from  a  crossflow  plane  (i.e.,  an  axial  station). 

•  Print  W,  'N'  planes  -  Specify  the  radial  and  circumferential  increments  for 
saving  the  ASCII  crossflow  plane  data  (i.e..  Write  out  the  data  at  every  ‘m’ 
radial  points  and  ‘n’  circumferential  points).  The  computational  effort  required 
for  J/0  can  be  minimized  by  not  writing  out  every  radial  and  circumferential  data 
point. 
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CHAPTER  5:  CODE  EXECUTION 


OVERVIEW 

This  section  is  used  to  select  which  cases  are  evaluated  as  well  as  what  type  of 
solution  is  performed.  A  step  by  step  procedure  for  executing  a  ran  is  provided  below 
along  with  a  sample  view  of  the  execution  screen. 

•  Select  the  desired  run. 

Note:  Multiple  runs  may  be  selected  in  a  number  of  ways.  The  first  is  to  hold 
down  the  'Ctrl'  key  and  click  on  each  of  the  desired  runs.  Alternatively, 
you  can  click  on  the  first  desired  run,  hold  down  the  ’shift'  key,  and  then 
click  on  the  last  run  to  select  all  cases  in-between.  Finally,  you  can  click  on 
the  ’select  AH'  button  to  run  all  available  cases. 

•  Select  the  type  of  ran  desired. 

Note:  If  'Full  Run'  is  selected  the  code  will  execute  the  Inviscid  Euler  solver, 
write  out  the  flowfield  at  the  axial  intervals  specified  in  the  Boundary  Laver 
Execution  Parameter  section,  and  then  execute  the  boundary  layer  solver  to 
get  an  approximate  viscous  solution.  If  'Inviscid  Run'  is  selected  the  code 
will  execute  the  Inviscid  Euler  solver  and  then  terminate.  No  viscous 
approximation  will  be  calculated.  If  a  viscous  calculation  is  desired,  click 
on  the  'Viscous  Run'  button,  and  select  the  desired  case.  Only  the  runs  that 
have  been  solved  with  the  inviscid  solver  will  show  up  as  choices  when  the 
'Viscous  Run'  button  is  selected.  This  is  because  the  viscous  boundary 
layer  solver  requires  the  inviscid  flowfield  as  a  boundary  condition. 

•  Specify  whether  a  symmetry  plane  exists. 

Note:  If  you  know  a  priori  that  the  geometry  is  non-symmetric  about  the  pitch 
plane,  you  must  select  'Full  Grid  (All  Runs)'.  If  the  geometry  is  symmetric 
about  the  pitch  plane  choose  either  ’symmetry  Plane  (AH  Runs)'  if  you  wiH 
not  be  running  any  asymmetric  cases,  or  alternatively,  choose  ’s)mimetry 
Plane  (Zero  RoH  and  Yaw)'  to  ran  a  symmetry  plane  only  for  those  cases 
with  zero  yaw  and  roH  angles.  If  a  symmetry  plane  is  being  considered  then 
only  half  of  the  computational  domain  (180  vs.  360  degrees)  is  considered. 

Running  with  a  symmetry  plane  will  cut  the  computational  time  approximately  in 

half. 

•  Click  on  'OK'  to  perform  the  desired  calculations. 
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•  After  the  runs  have  completed,  revisit  the  'Execution'  section  to  check  on  the 
termination  status  of  each  run.  If  the  'Runs'  column  contains  one  asterisk  (*)  then 
the  inviscid  solution  has  completed  (e.g.,  RunOOOl).  If  it  contains  double  asterisks 
(**)  then  both  the  inviscid  and  viscous  solutions  have  been  completed  (Run0002). 
In  addition,  the  'Exit  Status'  column  should  contain  'Normal  Exit'  if  no  runtime 
errors  were  encountered.  If  the  code  did  not  exit  normally  then  an  error  code  will 
be  provided.  These  error  codes  are  5000  series  numbers  for  die  Euler  solver,  6000 
series  numbers  for  the  boundary  layer  solver,  7000  series  numbers  for  the  cone 
solver,  and  8000  series  numbers  for  the  blunt  body  solver.  A  description  of  each 
of  the  possible  error  codes  is  presented  in  the  Error  Code  section. 


— . r 


RunOOOl 
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0.000 

0.000 

Run0002 

4.250 

0.000 

0.000 

Iw  1 « 

Run0003 

4.500 

0.000 

0.000 

Bjs  1  n 

RunOOO-^ 
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10.000 

0.000 
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RunOOOB 
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10.000 

0.000 
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RunOOOB 
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0.000 
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RunOOO? 

4.000 

20.000 

0.000 

0.000 
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20.000 

0.000 
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Run0009 
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20.000 

0.000 

0.000 

RunOOl 0 

4.000 
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0.000 

0.000 

RunOOIl 

4.250 

30.000 

0.000 

0.000 

RunOOl 2 

4.500 

30.000 

0.000 

0.000 

A 

Normal  Exit 

Me 

Normal  Exit 

Me 

Normal  Exit 

A 

Normal  Exit 

A 

Normal  Exit 

A 

Normal  Exit 

A 

5105 

HGURESl.  EXECUTIVE  CONTROL  PARAMETERS 

5.1  CONSOLE  OUTPUT 

5.1.1  Euler  Console  Output 

As  the  Euler  solver  is  executing,  a  variety  of  information  is  written  to  a  command 
prompt  (see  sample  below).  This  includes  messages  denoting  which  code  is  currently 
running  (cone,  blunt  body,  or  Euler  solver)  as  well  as  a  number  of  parameters  denoting  the 
status  of  the  run. 
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0 ;  \Users\df robins' 

\Zeus'\Distribution>echo  off 

7  Files's)  copied 

7  Files's)  copied 

Kxxx  (Run0001) 

xx  Free-Streatn  Starting  Solution 

>0 

xxxx  (Run0001)  xxs 

XX  Executing  ZEUS  (Free-Strea®  Inflow  Prof ilejixxxx 

STEP  =  25  ZETft 

-  .3.500E-02  DZETft  =  1.000E-03 

NCFL.HCFL  =  1 

18 

STEP  =  50  ZETft 

:  6.000E-02  DZETft  1.000E-03 

NCFL. MCFL  =  1 

G 

STEP  =  75  ZETft 

=  8,50OE-02  DZETft  =  1 .000E-O3 

NCFL.HCFL  =  1 

19 

NORHftL  STOP 

xxxx  (Run8001) 

XX  Run  Completed  xxxx 

Press  any  key  to 

continue  .  -  .  _ 

FIGURES!.  EULER  CONSOLE  OUTPUT 

The  output  written  to  the  command  prompt  for  the  inviscid  solver  includes  the 
following  parameters: 

•  Starting  profile  -  Denotes  which  type  of  starting  solution  is  employed 
(freestream,  cone,  or  blunt  body).  This  parameter  in  specified  in  the  Integration 
Control  section. 

•  Step  -  Counter  for  the  number  of  axial  steps  taken.  The  frequency  with  which 
the  screen  data  is  written  is  controlled  by  the  Output  Residual  parameter  in  the 
Output  Control  section. 

•  Zeta  -  Axial  location  of  the  current  integration  step. 

•  (Zeta  -  Local  axial  stepsize.  This  value  changes  as  the  solution  proceeds 
downstream  and  is  controlled  by  the  stability  criterion  specified  in  the 
Tntpfrfltinn  Control  section.  If  the  axial  stepsize  is  too  large  to  resolve  the 
geometric  features  then  die  maximum  allowable  stepsize  can  be  adjusted  in  the 
Litegration  control  section. 

•  nCFL,  mCFL  -  The  radial  and  circumferential  index  for  the  maximum  residual. 

•  Exit  Status  -  Upon  termination  the  exit  status  of  the  current  run  is  displayed 
along  with  any  error  messages  that  occurred  during  the  Euler  solution.  'Normal 
Stop'  is  displayed  if  no  runtime  errors  were  encountered. 

Note:  After  the  run  has  completed  the  window  will  display  'Press  Any 
Key'  to  continue.  When  a  key  is  depressed,  the  command  window  is 
deleted  and  control  is  returned  to  the  Zeus;*^  graphical  user  interface 
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window.  Whether  or  not  the  window  pauses  after  completion  is  controlled 
by  the  'Pause  After  Each  Run'  parameter  in  the  Options  section. 

5.1.2  Boundary  Layer  Console  Output 


s  C:\WINNT\System32\CMD.exe 

: \Users\df robins\Z9us\Distribytion>echo  off 
File(s)  copied 
?  File(s)  copied 


Run0801 )  KKKx 

Executing 

Boundary 

Lay^r 

Code  (ftdiabat 

Step  = 

108 

Zeta 

r 

6.SI860E-82 

Step  = 

280 

Zeta 

:: 

1 ,19?2E-01 

Step  = 

300 

Zeta 

r 

1  .7i|58E-01 

Step  = 

me 

Zeta 

r 

2,2Si1^E--81 

Step  = 

500 

Zeta 

:: 

2.8i|38E-01 

Step  = 

600 

Zeta 

:: 

3.3316E-01 

Step  = 

?06 

Zeta 

:: 

3.3^02E-01 

Step  = 

808 

Zeta 

= 

i|.4888E-01 

Step  = 

300 

Zeta 

5.03?^4E-01 

I  Normal  Stop 


FIGURE  33.  BOUNDARY  LAYER  CONSOLE  OUTPUT 

The  ouQiut  written  to  the  command  prompt  for  the  viscous  boundary  layer  solver 
includes  the  following  parameters: 

•  Case  description  -  During  each  run  the  case  number  being  executed  is  displayed 
(Runxxxx)  along  with  the  appropriate  heat  flux  boundary  condition  at  the 
surface  (i.e.,  Adiabatic  or  Isothermal).  The  boundary  condition  is  specified  in 
the  General  Aerodynamic  Data  section. 

•  Step  -  Counter  for  the  number  of  axial  steps  taken.  The  frequency  with  which 
the  screen  data  is  written  is  controlled  by  the  Output  Residual  parameter  in  the 
Output  Control  section. 

•  Zeta  -  Axial  location  of  the  current  integration  step. 

•  Exit  Status  -  Upon  termination  the  exit  status  of  the  current  run  is  displayed 
along  with  any  error  messages  that  occurred  during  the  Euler  solution.  'Normal 
Stop'  is  displayed  if  no  runtime  errors  were  encountered. 

Note:  After  the  run  has  completed  the  window  will  display  'Press  Any 
Key'  to  continue.  When  a  key  is  depressed,  the  command  window 
is  deleted  and  control  is  returned  to  the  Zeus"^  graphical  user 
interface  window.  Whether  or  not  the  window  pauses  after 
completion  is  controlled  by  the  'Pause  After  Each  Run'  parameter  in 
the  Options  section. 
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5.2  ERROR  MESSAGES 
Overview 

The  possible  error  messages  encountered  while  executing  the  Zeus"*^  code  are  hsted 
below  for  each  section  of  the  tool.  Note  that  5000  series  errors  occur  in  the  Euler  solver, 

6000  series  in  the  boundary  layer  solver,  7000  series  in  the  cone  solver,  and  8000  series  in 
the  blunt  body  solver.  A  description  for  each  series  is  listed  in  the  corresponding  section 
below. 

Zeus  Code 

Rnnndary  Laver  Code 
Cone  Code 
Blunt  Body  Code 

5.2.1  Zeus  Error  Messages 

Maximum  iteration  count  to  locate  comer  #l-#4  exceeded. 

Negative  pressure  encountered  while  decoding  conservative  variables. 
Undefined  surface  geometry  parameter. 

Forward  swept  leading  edges  (control  surfaces)  not  yet  implemented. 
Rearward  swept  trailing  edges  (control  surfaces)  not  yet  implemented. 

After  first  order  fix  the  flow  remains  subsonic. 

Negative  square  root  occurred  in  the  Prandtl-Mayer  expansion  calculation. 
Error  in  linear  interpolation  for  shockwave. 

Dimensions  are  too  large. 

MAZ()  array  not  consistent  with  MA()  array. 

Specified  radial  clustering  not  attainable  for  aft  GASP  grid. 

Undefined  control  surface  input. 

Clustering  parameters  not  consistent  between  the  Cone  code  and  the  Zeus 
code. 

Error  in  oblique  shock  calculation. 

No  possible  oblique  shock  solution. 

No  convergence  in  Prandtl-Meyer  calculation. 

The  number  of  specified  axial  plotting  stations  is  greater  than  the  maximum 
allowed.  The  MaxPlots  parameter  must  be  increased  in  the  parameter 
statement  and  the  code  must  be  recompiled. 

Invalid  starting  solution  passed  to  the  Euler  solver.  Typically  caused  by 
errors  in  the  cone  or  blunt  body  starting  profile  generators. 

5.2.2  Boundary  Laver  Error  Messages 

6001  Maximum  iteration  count  reached 

6002  Invalid  axial  force  calculation 


5101-5104 

5105-5106 

5108 

sm 

5110 

5111-5112 

5115-5116 

5118 

5120 

5121 

5123 

5124 

.512R-5129 

51.30 

5131 

5132 
5134 


5135 


Error  messages  for  the  Zeus  Euler  code  (inviscid 
solution). 

Error  messages  for  the  Boundary  Layer  code  (viscous 
solution). 

Error  messages  for  the  Cone  code  (generates  starting 
profiles). 

Error  messages  for  the  Blunt  Body  code  (generates 
starting  profiles). 


49 


NSWCDD/rR-98/147 


5.2.3  Cone  Error  Messages 

7001  Shock  angle  too  great.  If  excessive  turning  angles  are  encountered  during 

the  cone  solution  then  the  initial  profile  will  be  subsonic. 

5.2.4  Blunt  Body  Error  Messages 

8001  Errors  in  the  blunt  body  starting  profile  generator. 
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CHAPTER  6;  POST-PROCESSOR 


OVERVIEW 

The  Zeus‘S  tool  includes  a  built-in  post-processor  for  viewing  the  solution  data.  A 
variety  of  methods  of  examining  the  results  are  available  including  ASCII  representations 
of  the  data  as  well  as  plots  of  a  wide  variety  of  parameters,  all  of  which  are  discussed  in 
detail  below.  The  plotting  features  of  the  Zeus^  code  make  use  of  the  Tecplol^  plotting 
package. 

Note:  Zeus"^  can  be  used  if  Tecplot  is  not  installed  on  the  users'  machine, 
however,  the  plotting  capabilities  of  the  post-processor  will  not  be 
functional.  ASCII  representations  of  the  data  will  still  be  available. 


The  'Postprocessor'  menu  option  offers  the  following  commands: 


ASCII  Forces/Moments 


Single  Run  Plots 


Multi-Run  Plots 


Examine  an  ASCII  data  file  of  the 
forces/moments  and  centers  of  pressure  acting  on 
the  body. 

T  .inp.  plots  of  the  surface  pressures,  forces,  and 
moments.  Vector  and  streamline  plots.  Contour 
plots  of  a  variety  of  flowfield  parameters  (Mach 
number,  density,  temperature,  ...etc.). 
Comparative  line  plots  of  forces/moments  versus 
a  variety  of  independent  variables  (Mach  number, 
angle  of  attack,  yaw  angle,  and  roll  orientation). 


6.1  ASCn  Data 
6  .11  Forces  &  Moments 

There  are  two  options  when  examining  the  forces,  moments,  and  centers  of 
pressure  for  a  given  set  of  runs.  The  first  is  to  view  the  totals  or  the  sum  from  each  of  the 
zone  edges.  This  gives  the  sum  of  the  values  from  the  missile  surface  and  the  control 
surfaces,  or  the  total  acting  on  the  geometry.  The  second  option  prints  the  forces/moments 
separately  on  each  edge  of  each  zone.  This  is  often  useful  if  the  values  on  a  particular 
control  surface  are  desired  (i.e.,  force  due  to  1  fin).  Recall,  that  edge,  refers  to  the  mssile 
surface,  edgej  and  edgej  refer  to  the  zonal  boundaries  (i.e.,  the  control  surfaces,  if  they 
exist),  and  edge4  refers  to  the  outer  boundary  (see  figure  below). 
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FIGURE  34.  EULER  GRID  REFERENCE  DIAGRAM 

For  each  case  run,  the  ASCII  output  file,  'ForcesMoments.txt’  will  contain  the 
following  data: 

Mach# 

Angle  of  Attack 
Angle  of  Yaw 
Roll  Angle 

Normal,  Side,  and  Axial  Force 
Yaw,  Pitch,  and  Roll  Moments 
Xcp/Ljjgf  and  Y^p/L^gf  (centers  of  pressure) 


Note:  An  approximate  value  for  the  axial  force  due  to  the  base  (i.e.,  base  drag)  is 
calculated  using  the  methodology  of  reference  14.  Note  that  the  corrections 
due  to  angle  of  attack  and  fin  thickness  have  been  omitted.  Also,  note  that 
the  base  drag  calculation  requires  the  ratio  of  the  base  area  to  the  reference 
area  Zeus‘S  calculates  the  base  area  assuming  a  circular  cross- 

section,  and  therefore,  if  a  non-circular  cross-sectional  geometry  is  applied, 
the  resultant  base  drag  must  be  multiplied  by  the  factor  .  In  other 

words,  an  equivalent  circular  base  area  must  be  applied. 

If  both  the  inviscid  and  viscous  solutions  have  been  performed  then  the  ouqrut  for 
the  forces/moments  will  be  listed  componentwise.  In  other  words,  the  inviscid,  viscous, 
and  induced  contributions  will  be  listed  along  with  the  totals. 
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6.2  SINGLE  RUN  PLOTS 
Overview 

This  section  of  the  post-processor  is  used  to  examine  the  output  data  for  a  single 

run.  To  make  a  comparative  analysis  of  different  runs  see  the  M^ti^ffllHots  section. 

The  ’single  Run  Plots'  menu  offers  the  following  three  choices: 

Surface  Pressures  Line  plot  of  the  surface  pressures. 

Forces  and  Moments  Line  plots  of  the  forces,  moments,  and  centers  of 

pressure. 

ContourA^ector  Plots  Vector  and  streamline  plots.  Contour  plots  of  a 

variety  of  flowfield  parameters  (Mach  number, 
density,  temperature,  ...etc.). 

6.2.1  Surface  Pressures 

The  first  option,  'Plot  Surface  Pressure',  will  display  a  graph  of  the  surface 
pressure  versus  axial  distance  for  zone  edges  1,  2,  and  4  (see  figure  in  A$Cn_d&ta 
section).  The  pressures  are  listed  for  each  circumferential  plane  on  the  missile  surface  and 
for  each  radial  plane  on  the  zone  edges  (i.e.,  control  surfaces).  Each  plane  is  labeled  as 
'Plane  #x.  Zone  #y.  Edge  #z'.  If  the  Edge  number  is  '  1 '  then  the  Plane  #  refers  to  the 
circumferential  index.  If  the  Edge  number  is  '2'  or  '4'  then  the  Plane  #  refers  to  the  radial 
index. 

6.2.2  Forces  and  Moments 

The  second  single-run  plotting  option,  'Plot  Forces  &  Moments',  will  plot  either 
the  total  force/moment  or  the  edge  force/moment  as  a  function  of  axial  location.  The  totals 
versus  edge  values  are  described  in  the  ASCII  data  section.  The  'x',  'y ',  'z'  coordinates 
refer  to  normal,  side,  and  axial  components,  respectively. 

6.2.3  ContourA^ector  Plots 

This  plotting  option  generates  contour  plots  of  a  variety  of  flowfield  variables  for  a 
given  run.  Click  on  'Choose  Plot  Attributes'  to  bring  up  the  dialog  box  shown  below  and 
specify  the  plotting  attributes.  Select  the  dimensions  and  the  variables  desired  for  the 
contour  plot.  Note  that  the  coordinates  (x,  y,  z)  and  the  velocities  (u,  v,  w)  are  always 
written  to  the  contour  plot  file.  After  the  plot  information  has  been  specified,  click  'OK'  to 
accept  the  values.  Then  select  'Plot  Data'  to  generate  the  contour  plots. 


53 


NSWCDD/TR-98/147 


HGURE  35.  CONTOUR  PLOT  CONTROLS 


6.3  MULTI-RUN  PLOTS 
Overview 


This  section  of  the  post-processor  is  used  to  compare  the  effect  of  Mach  number, 
angle  of  attack,  yaw  angle,  and  roll  orientation  on  the  forces,  moments,  and  centers  of 
pressure  of  a  given  geometry.  Select  the  desired  independent  variable  as  well  as  the  desired 
dependent  variables  (see  figure  below).  Click  on  'OK'  to  generate  a  plot  for  all  of  the 
completed  runs. 
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nGURE36.  24  PLOT  CONTROLS 
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A.1  INTRODUCTION 


STARTING  ZEUS++ 

Execution  of  the  code  begins  by  double  clicking  on  the  program  'ZeusPP.exe’,  entering  a 
name  for  the  current  case  ('tutorial'),  and  then  clicking  on  ’start  New  Case'.  For  the  tutorial,  we 
will  be  considering  a  SeaSparrow  missile^  '^  with  both  dorsal  and  tailEns. 

Note:  The  name  entered  cannot  contain  blank  spaces. 

Execution  of  the  Zeus*’^  code  can  be  divided  into  four  basic  sections: 

1.  Geometry  Setup 

2.  Case  Inputs 

3.  Code  Execution 

4.  Post-Processing 


A.2  GEOMETRY  SETUP 


OVERVIEW 

The  geometry  setup  section  can  be  sub-divided  into  four  areas: 

1)  Specify  the  grid  parameters 

2)  Define  the  missile  geometry 

3)  Define  the  control  surfaces 

4)  Generate  the  computational  domain  (optional). 


A.2. 1  GRID  PARAMETERS 

The  properties  that  define  the  computational  grid  are  specified  in  this  section.  Each  of 
these  parameters  is  described  in  detail  in  the  Grid  Parameters  section. 
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•  From  the  main  menu,  click  on  'Geometry'  and  then  select  'Grid  Parameters'  and  'Euler 
Grid  Properties.  For  an  explanation  of  each  of  the  parameters  on  this  screen,  press  the 
'FI'  key.  The  default  values  are  adequate  for  the  current  tutorial,  therefore,  simply 
click  on  'OK'  to  accept  the  current  values. 

Note:  From  the  'Grid  Parameters'  option,  it  is  seen  that  there  are  two  sets  of  grid 
parameters  available.  In  the  first  (Euler  Grid  Properties),  the  parameters  specified 
are  applied  during  both  the  Zeus  Euler  solution  as  well  as  the  viscous  boundary 
layer  solution.  In  the  second  (Navier-Stokes  Grid  Properties),  the  parameters 
specified  are  used  to  generate  a  grid,  which  can  be  exported  to  a  full  Navier-Stokes 
computational  fluid  dynamics  solver.  The  Navier-Stokes  parameters  are  used  only 
for  grid  generation  purposes,  and  therefore,  if  the  user  does  not  intend  to  export 
the  grid  to  an  N-S  solver,  this  section  may  be  omitted. 


A.2.2  MISSILE  SECTIONS 

After  specifying  the  grid  parameters,  the  next  step  is  to  define  the  missile  geometry.  For 
the  tutorial,  the  missile  geometry  will  be  entered  as  three  separate  sections:  1)  Tangent-Ogive 
nose,  2)  cylindrical  mid-body  and  3)  boat-tail  flare. 

The  first  section  is  the  Tangent-Ogive  Nose 

•  Click  on  'Geometry'  and  then  select  'Missile  Sections'. 

•  Click  on  the  down  arrow  next  to  'select  Geometry  Type'  and  select  'Tangent  Ogive'. 

•  Click  on  'Add  New  Section'. 

•  Enter  the  final  axial  location  (0.0686)  of  the  nose  as  well  as  the  final  radius  (.01524). 

Note  that  the  missile  starts  at  0,0. 

•  Click  on  'OK'  and  the  code  will  calculate  and  display  the  initial  missile  section. 

The  next  section  is  the  cylindrical  mid-body  and  is  described  by  a  conic  section. 

•  Click  on  the  down  arrow  next  to  'select  Geometry  Type'  and  then  select  'Conic 
Section'. 

•  Click  on  'Add  New  Section'. 

•  The  initial  axial  location  and  radius  are  copied  from  the  end  of  the  previous  section. 

•  Enter  the  final  axial  location  (0.537)  of  the  conic  as  well  as  the  final  radius  (.01524). 

•  Click  on  'OK'. 

The  final  section  is  the  boat-tail  flare  and  is  described  by  a  conic  section. 

•  Click  on  the  down  arrow  next  to  'select  Geometry  Type'  and  then  select  'Conic 
Section'. 
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•  Click  on  'Add  New  Section'. 

•  The  initial  axial  location  and  radius  are  copied  from  the  end  of  the  previous  section. 

•  Enter  the  final  axial  location  (0.5486)  of  the  conic  as  well  as  the  final  radius 
(.013715). 

•  Click  on 'OK'. 

•  Click  on 'Finished' 

•  Press  'Ctrl-R'  to  regenerate  the  geometry. 


A.2.3  CONTROL  SURFACES 

After  specifying  the  missile  geometry,  the  next  step  is  to  define  the  control  surfaces. 
For  the  tutorial,  two  sets  of  control  surfaces  are  considered.  The  first  is  the  dorsal  fins  and  the 
second  is  the  tail  fins.  The  control  surfaces  must  be  entered  in  the  order  in  which  they  appear 
on  the  missile,  therefore,  the  dorsal  fins  must  be  entered  first,  followed  by  the  tail  fins. 

The  first  set  of  control  surfaces  (dorsal  fins)  will  be  generated  using  the  'Modified 
Wedge'  option. 

•  Click  on  'Geometry'  and  then  select  'Control  Surfaces'. 

•  Click  on  the  arrow  next  to  ’select  Control  Surface'  and  select  'Modified  Wedge'. 

•  Click  on  'Add  Control  Surface'. 

•  Enter  a  label  for  the  control  surface  (“Dorsal  Fins”). 

•  Enter  the  axial  location  of  the  control  surface  hinge  line  (“0.2743”).  The  hinge  line  is 

the  point  on  the  control  surface  about  which  deflections  are  made. 

•  Enter  the  offset  angles  (in  degrees)  for  each  of  the  fins.  For  the  current  case  the  fins 

are  located  in  the  '+'  configuration,  and  therefore,  the  offset  angles  are  90,  180,  270, 
and  360°. 

Note:  It  is  possible  to  delete  an  individual  fin  by  right  clicking  on  the  fin  you  wish  to 

remove.  Both  the  offset  angle  and  the  deflection  angle  will  be  replaced  with  ' - '  to 

indicate  that  they  have  been  eliminated.  In  order  to  add  the  fin  back  into  the 
calculation,  simply  right  click  on  it  a  second  time.  This  feature  is  useful  when  you  want 
to  have  a  different  number  of  fins  for  each  set  of  control  surfaces  (i.e.,  2  dorsal  fins  and 
4  tail  fins). 

•  Enter  the  deflection  angles  (in  degrees)  for  each  of  the  fins.  Use  +5.0  for  the  90  and 

270°  fins. 

Note:  The  fins  are  always  deflected  in  a  clockwise  orientation  (i.e.,  entering  +5  for  both 
fins  will  generate  an  asynunetric  deflection). 

•  Click  'OK' 
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•  The  parameters  used  to  describe  the  modified  wedge  control  surface  are  described  in 
detail  in  the  Modified  Wedge  section.  For  the  Dorsal  fins,  enter  the  following  values 
and  click  on  'OK'  when  completed: 


1  ’  Modified  Wedge 

Control  Surface 

HIslE3| 

Enterthe  leading  edge  half  angle 

[JoTs 

::[deg] 

Enterthe  trailing  edge  half  angle 

|3.977 

[deg] 

\  Enterthe  leading  edge  sweep  angle 

Hs 

:[deg] 

Enterthe  span  (rootto  tip  distance) 

J'0597 

Distance  from  leading  edge  to  hinge  line 

[.0422 

Root 

Tip 

Enterthe  leading  wedge  thickness 

[.002225 

;  |3.822e-4 

Enter  the  trailing  wedge  thickness 

(.002225 

;  |5.04e-4 

Enterthe  chord 

[0742 

[.0145 

OK 

Cancel  | 

Convert  to ‘User  Specified  Geometiy*  | 

X' 


FIGURE  A- 1.  MODIFIED  WEDGE  SAMPLE  INPUT 


The  tailfms  represent  the  second  set  of  control  surfaces.  Because  of  the  layout  of  the 
geometry  of  these  fins,  they  will  be  generated  using  the  'User  Defined  Geometry'  option. 


•  Click  on  the  arrow  next  to  ’select  Control  Surface'  and  select  'User  Defined  Geometry'. 

•  Click  on  'Add  New  Section'. 

•  Enter  a  label  for  the  control  surface  (“Tail  Fins”). 

•  Enter  the  axial  location  of  the  control  surface  hinge  line  (“0.5167”). 

•  Enter  the  offset  angles  (in  degrees)  for  each  of  the  fins.  For  the  current  case  the  fins 
are  located  in  the  ’+'  configuration,  and  therefore,  the  offset  angles  are  90,  180,  270, 
and  360°. 

•  Enter  the  deflection  angles  (in  degrees)  for  each  fin  (0°) 

•  Click  'OK' 

•  To  ensure  that  the  tail  fins  are  described  by  planar  surfaces  the  geometry  must  be 
entered  as  three  cross  sectional  cuts.  A  detailed  overview  of  the  methodology  used  to 
generate  this  control  surface  geometry  can  be  found  at  User  Defined  Geometry.  For 
the  tail  fins  under  consideration,  enter  the  following  values  for  the  three  cross- 
sectional  planes: 
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BSQI 


(Plane#!  of 3) 


Enteracross-seclionaj  plane  of  the  ;  'X*  points 

/:  cxintrol  surface.  These  points  are  given  ■  -v  ^  ' '  " '  JZ 

f  .with  respectto  the  hinge  line  (i.e.z-0  is  Vpoints  lO.OOOOO  0.00000  0.00000 .001 525 
I  the  axial  location  of  the  hinge  line)  *  ^  - 


7’  points 


|- 0390 -.01 86545  .011 5545  .031 9 

(Qse  spaibeS  between  points) 


Next  Plane 


Delete  Rane 


Rnished 


FIGURE  A-2.  USER  DEFINED  CONTROL  SURFACE  SAMPLE  INPUT  (PLANE  1) 


Click 'Next  Plane' 


:  Control  Surface  Definition 


(Rane  #2  of  3) 


k  Enter  a  cross-sectional  plane  of  the 
I  coiitrol  surface.  These  pointe  are  given 
J  with  respect  to  the  hinge  line  (i.e.  z»0  is 
S;  the  axial  location  of  the  hinge  line) 


’X'  points 


Vpoints  |0.01 968  001^ 0.01 968 .021 205 

•2.  points  |- 008791  .0115545 .0115545 .0319 

(Use  spaces  between  ppints) 


Next  Rane 


Delete  Rane 


Rnished 


FIGURE  A-3.  USER  DEFINED  CONTROL  SURFACE  SAMPLE  INPUT  (PLANE  2) 


Click  'Next  Plane' 
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:  Controf  Surface  Definition 


mmm 


(Rane#3  of  3) 


0.0  0  0  0.0  0.0 


Enter  a  cross-sectional  plane  of  the  ‘X' points  ’ 
y  control  surface.  These  points  are  given  ~  r— — — — — rr — 

J  with  respect  to  the  hinge  line  O.e.  z»0  is  y‘  points  [  o^7725  0-^7725.0-17725  04  7725 
f  the  axial  location  of  the  hinge  line)  * 


‘2’ points  |.0319.0319.0319.0319 

(Lise  spaces  between  points) 


Next  Plane 


Delete  Rane 


Rnished 


FIGURE  A-4.  USER  DEFINED  CONTROL  SURFACE  SAMPLE  INPUT  (PLANE  3) 

•  Click  'Finished' 

•  Click  'Finished' 

•  Press  'Ctrl-R'  to  regenerate  the  geometry 


A.2.4  GRID  GENERATION  (OPTIONAL) 

•  View  Euler  Grid 

The  purpose  of  this  section  is  to,  examine  the  inviscid  grid  that  will  be  used  by  the  Zeus 
Euler  solver  to  ensure  that  an  acceptable  grid  is  being  employed. 

Note:  1)  The  axial  step  size  for  the  inviscid  solution  is  calculated  at  runtime,  however, 
the  axial  step  size  for  grid  generation  purposes  is  taken  to  be  constant. 

2)  The  outer  boundary,  for  grid  generation  purposes,  is  taken  to  be  that  specified 
in  the  Euler  grid  parameter  section. 

•  Click  on  'Geometry'  and  then  select  'Grid  Generation'  and  'View  Euler  Grid' 

•  Click  'Plot  Grid  as  a  Series  of  2d  planes' 

•  Enter  a  starting  axial  location  (0.01) 

•  Enter  a  final  axial  location  (0.5486) 

•  Enter  the  number  of  axial  step  (25) 

•  Click  'OK' 

•  After  the  Tecplot  window  appears: 

•  Click  'File',  'Macro',  and  then  'Play' 

•  Double  click  on  the  ’2dGrid_Movie.mcr'  file  in  the  'Playback'  directory. 
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Note:  The  macro  file  will  step  axially  down  the  missile  and  take  a  snapshot  of  each 
location.  The  viewport  shown  during  the  movie  generation  (i.e.,  x  and  y-axis 
ranges)  can  be  specified  in  Zeus**  before  generating  the  Euler  grid  by  clicking  on 
'View'  and  then  selecting  'Options'.  Enter  the  desired  values  for  the  'x'  and  'y' 
axis  ranges.  A  detailed  description  of  all  of  the  parameters  is  provided  in  the 
View/Options  section. 

•  A  movie  of  the  grid  will  be  created  and  written  out  as  'Playback\Zeus_grid2d.rm . 

This  movie  file  can  be  replayed  using  the  Tecplot  Framer  program. 

•  Navier-Stokes  Grid 

The  Zeus**  code  functions  both  as  an  inviscid  solver  (with  approximate  viscous 
corrections),  and  a  CFD  grid  generator.  This  allows  one  to  export  a  three-dimensional 
computational  grid  and  the  appropriate  boundary  conditions  to  a  full  Navier-Stokes  solver. 
Currently,  the  grid  and  boundary  condition  files  generated  are  suitable  for  the  GASP  flow 
solver.  One  circumferential  plane  of  the  Seasparrow  Navier-Stokes  grid  is  shown  below  with 
an  explanation  of  each  parameter  provided  in  the  Navier-Stokes  Grid  section. 


FIGURE  A-5.  SAMPLE  OUTPUT  OF  NAVIER-STOKES  GRID  GENERATOR 
(2d  CROSS-SECTION  SHOWN) 

•  Click  on  'Geometry',  'Grid  Generation',  'Navier-Stokes  Grid',  and  'Generate  Grid' 

•  Enter  the  number  of  axial  sections  (6).  The  geometry  is  separated  into  a  number  of 
axial  stations  so  that  clustering  can  be  performed  at  a  number  of  axial  locations.  From 
the  above  figure,  the  six  axial  sections  are  shown  between  the  vertical  lines.  The  six 
sections  are  as  follows: 

1 .  Nose  section 

2.  Section  between  nose  and  first  set  of  fins 

3.  Section  covering  dorsal  fins 

4.  Section  between  dorsal  fins  and  tail  fins 

5.  Section  covering  the  tail  fins 

6.  Section  aft  of  the  missile 
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Note:  If  you  extend  the  grid  beyond  the  aft  end  of  the  missile  (shown  in  blue  between 
the  vertical  lines),  then  an  additional  section  will  be  included  to  calculate  the 
flow  directly  behind  the  missile  (shown  in  yellow  between  the  horizontal  lines). 

•  Enter  the  number  of  axial  stations,  ending  axial  location,  beginning  spacing,  and 
ending  spacing  for  each  of  the  six  sections  (values  provided  below): 

Note:  To  simplify  the  initial  estimate  for  axial  spacing  a  'Reset  Spacing'  button  is 
provided.  Simply  enter  the  number  of  axial  stations  and  the  ending  axial 
location  for  each  of  the  desired  sections.  Then,  click  on  'Reset  Spacing'  to  obtain 
an  equally  spaced  axial  grid.  This  grid  can  then  be  modified  to  provide  the 
desired  axial  clustering. 


:  Navier-Stokes  Grid  Properties 


-Axial  Clustering 

Enter  the  number  axial;  sedidns  W 


i  Reset  Spacing 


#of  T  stations 


Ending  'Z'  location  Beginning  Spacing  Ending  Spadng 


20  ;  ; 

.0686 

.003 

12  :■ 

.230 

.01  If, 

20  ;; 

.3064 

.003  ■ 

12 

.4775 

.003  . 

21  i  ^  . 

.5486  •  > 

.003 

20  r:  . 

.7 

.001  ;  ;;r: 

Aft  Grid  - 

P  Include  grid  aft  of  missile 

Enterthe  #  of  radial  pts.  fisf 

(centerline  to  surface) 


-Outer  Boundaiy  Clustering  — 
j  Centerline  to  first  point  angle 
Centerline  to  last  point  angle  lii” 


Cancel 


FIGURE  A-6,  SAMPLE  INPUTS  FOR  NAVIER-STOKES  GRID  GENERATOR 


•  For  the  aft  grid  section,  select  the  'Include  Grid  Aft  of  Missile'  and  enter  the  number  of 
radial  points  to  include  (15). 

•  In  order  to  control  the  clustering  of  the  outer  boundary  of  the  nose,  the  beginning  and 
ending  angles  (see  'a'  and  'b'  in  the  above  figure)  of  the  outer  boundary  are  specified. 
Enter  the  values  shown  in  the  above  figure. 

•  Click  on  'Ok'  to  bring  up  the  boundary  condition  window  (see  Figure  A-7). 

•  Select  the  surface  boundary  condition. 

•  Select  the  two-equation  turbulence  model. 

•  Click  on  'OK'  to  generate  the  computational  grid  and  boundary  condition  files. 
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Note:  The  output  from  the  grid  generation  code  is  placed  in  the  'filename/RunGrid' 
directory.  For  the  current  case,  this  is  'Tutorial/RunGrid'.  The  two  files  are  the 
grid  geometry  in  PlotSd  format  (GASPGrid.pSda)  and  the  GASP  boundary 
condition  input  file  (GASPbc.inp). 


■  Navier-Stokes  Boundaiy  Conditions 


^Surface  Boundaiy  Condition - 
Invisdd  No  Slip -Adiabatic  P^No  Slip -Isothermal 


p Two-Equation T urbulent Model  v>--- — 

Note:  The  type  of  two-equation  modal  applied  during  the  NS  solution 

must  be  specified  in  the  GASP  boundaiy  Gonditioii  input  file. 

(*  ^o  Two-Equation  Turijulence^^M^^ 

C  High  Reynolds  number  K-Epsilon 

r  LafTvBremhorsts  low  Reynolds  number  KrEpsilon  model 
P  Ghien*8  low  Reynolds  number  K-Epsilon  mode! 

P  \Wlcox's  KrOmega  model 

P  .Wilcox’s  low^eynolds  number  K-Omega  mode! 


i  , 

I  ' 

M 


Generate  Grid 


Generate  &  View 


Cancel 


FIGURE  A-7.  SAMPLE  BOUNDARY  CONDITION  INPUTS  FOR 
NAVIER-STOKES  GRID  GENERATOR 


A.3  INPUT  PARAMETERS 


OVERVIEW 

After  the  computational  doniain,  missile,  and  control  surface  geometries  have  been 
specified,  the  next  step  is  to  input  the  following  parameters: 

1.  Aerodynamic  Data 

2.  Integration  Controls 

3.  Separation  Modeling  Parameters 

4.  Boundary  Laver  Controls 

5.  Output  Controls 
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A.3.1  AERODYNAMIC  DATA 

This  section  is  used  to  specify  the  Mach  numbers  and  missile  orientations  (angle  of 
attack,  yaw,  and  roll)  to  be  used  during  execution  of  the  flow  solver.  Rather  than  specifying 
singular  values,  sweeps  are  defined  for  each  of  the  above  mentioned  parameters.  In  order  to 
generate  a  given  sweep,  input  the  beginning,  ending,  and  incremental  value  for  the  desired 
parameter  (see  below). 

•  Click  on  'Inputs'  and  then  select  'Aerodynamic  data'. 

•  Select  the  'Aerodynamic  Run  Matrix'  tab  and  then  enter  the  values  shown  in  the  figure 
below. 


I  Aerodynamic  Run  Matrix  )  General  Conditions  1 

|:  .  . . v:;' 


-Mach  Number  Sweep  — 

p  Configuration  Roll  Position  Sweep  -  -  ^ 

Rrst  Mach  Number 

■ 

Rrst  Roll  Position 

|o  1 

■ 

Last  Mach  Number 

I;- 

Last  Roll  Position 

\ 

Mach  Numbe  r  Interval  s 

Roll  Position  Intervals 

—  Annin  nf  AftAnL*  ^ 

■ft 

rviyic  Ul  wWbc?P 

Rrst  Angle  of  Attack 

LieTiecuon  oweep 

.  First  Tail  pefiection 

.s  ■ 

Last  Angle  of  Attack 

|20  .  .. 

L&sf  Tail  Dsfl action 

-J  • 

Angle  of  Attack  Inteivals 

10  I; 

^  Tail  Dati action  Intervafs 

. . .  ■ 

•ii 

-Angle  of  Yaw  Sweep — -:■■■  -  :  - — — 

f  - 

Rrst  Angle  of  Yaw  0  F 

LastAngledfYaw  \4  p 

Default  Values 

i; 

Angle  of  Yaw  Intervals  A 

OK 

;  Cancel 

"  Apply  ; 

FIGURE  A-8.  SAMPLE  OF  AERODYNAMIC  RUN  MATRIX  PARAMETERS 

The  next  section  is  used  to  specify  the  remaining  aerodynamic  parameters  including 
free-stream  atmospheric  conditions,  units,  reference  conditions  (diameter,  area,  and  moment 
arm),  viscous  parameters  (wall  heat  flux  boundary  condition,  Prandtl  number),  as  well  as  the 
gas  type  (real  or  ideal  gas). 

•  Click  on  the  'General  Conditions'  tab  to  specify  the  remaining  aerodynamic 
parameters. 
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•  The  free-stream  aerodynamic  parameters  (density,  pressure,  and  ratio  of  specific  heats) 
can  be  specified  in  a  variety  of  ways.  For  the  current  example,  a  unit  Reynolds 
number  is  provided  and  the  free-stream  quantities  are  computed.  Alternate  methods 
include  explicitly  specifying  the  values  or  specifying  an  altitude  and  using  standard 
atmospheric  tables  to  look  up  the  corresponding  values. 

•  Specify  the  units  as  metric. 

•  Enter  the  reference  diameter,  area,  and  moment  arm  shown  in  the  figure  below. 

•  Select  'Adiabatic'  and  specify  a  Prandtl  number  Of  0.72. 

Note:  The  viscous  parameters  are  not  used  in  the  Zeus  Euler  solver.  These  values  are, 
however,  required  inputs  for  the  approximate  viscous  (boundary  layer)  solver. 
Therefore,  these  values  are  not  required  unless  you  are  running  the  viscous 
solution  along  with  the  inviscid  flow  solver. 

•  Click  on  'OK'  and  wait  for  the  code  to  generate  the  required  run  matrix. 


f  i^rodyriamic  F^n  Mabw  General  Conditions 


- Free^ecim  Data — — 
r.{I962;U.S.  Standard  Tables 
rT959V^CT«ibles 

y^titude  10  i  [ft 


C,  fissure  |8T1 7  27 

Densit/  071 


- UnKs."  ^ .  ^ 

r  .English  (p-lbf/rt  rho-slug/fra 

Metric 


Gamma  |l 


<5*  Unit  Reynolds  Number 
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A.3.2  INTEGRATION  CONTROL 

This  section  is  used  to  specify  integration  parameters,  all  of  which  are  described  in 
detail  in  the  Integration  Parameters  section. 


’  Click  on  'Inputs'  and  then  select  'Integration  Control' 

'  Enter  all  of  the  values  shown  in  the  figure  below  and  click  'OK'. 


'  Integration  Parameters 


p  Space  Marching 


HSB 


Adal  location  to  begin  computation  |0.01 

Mel  location  to  terminate  computation. 

Step  number  to  terminate  corhputation. 


0.5^86 


25000 


-Starting  Solution  — 
r  Free  Stream 
^  f^ne  SQlutibr^ 
r-BluntBody 
R  Outer  boundary  tracks  shock 


Set  ParaiTieterB 


p.CFL------ - - - - - 

Step  size  safe^  factor  [0.<FCFL<1 .] 
;;PypiGallyuse0.9] 

■JO.9  1; 

Maximum  allowable  step  size 

■  .001 

-Riemann  Solver- 


C  Complete  Solver 

Approximate  Solver 
[Must  use  for  real  gas] 


-Bow  Shock  Flux- 


(5*  Use  Complete  Riemann  Solver 
:  t  lSHFLUX  -  0] 

C  Use  Freestream  Properties 
[ISHFLUX  -  1]^^-  ^ 


pBowShockAngle- 


r?  Use  Complete  Rierhann  Solver 
[ISHAMG  »^0] 

C  Use  Roe  Averaged  Variables 
[iSHANG-l] 


p  Interior  Point  UmltBr— - — - — — 

Umiterfor  points  interior  to  zones  [0.<XKI<2.]. 


p Turning  Angle  Correction 
Gorrectipn  factor 
[P<dfac<1.0] 


0.9 


OK 


Default  Values 


Cancel 


FIGURE  A-10.  SAMPLE  INPUT  FOR  INTEGRATION  PARAMETERS 


A.3.3  SEPARATION  MODELING 

This  section  is  used  to  specify  separation-modeling  parameters,  all  of  which  are 
described  in  detail  in  the  Separation  Modeling  section.  Separation  modeling  is  not  used  for  the 
tutorial,  and  therefore,  this  section  may  be  omitted. 


A.3 .4  BOUNDARY  LAYER  DATA 

If  an  approximate  viscous  solution  is  desired  then  the  following  values  must  be 
specified.  An  approximate  viscous  solution  is  obtained  by  solving  the  integral  form  of  the 
boundary  layer  equations.  The  integral  boundary  layer  solution  employed  requires  the  in  viscid 
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solution  from  the  Zeus  Euler  solver  as  a  boundary  condition.  An  in-depth  description  of  all  of 
the  parameters  shown  below  are  described  in  the  Boundary  Laver  Parameters  section. 

•  Click  on  'Inputs'  and  then  select  'Boundary  Layer  Data'. 

•  Enter  the  values  shown  in  the  figure  below  and  click  on  the  'Modeling  Parameters'  tab. 


I  Execution Porometere  | Modeling PareLmeters]  Ou^ut Parameters] 
i  pBoundaiy  L^r  Parameters  ^ 

Output  iiTviscldflowfield  every ‘iplotbroxiol  steps  1000 
Output  inviscid  flowfield  if  axial  stepsize  >  'delbl'  999.9 


Integration  Control — ~ — 
Starting  axial  location  for  integration 

Minimum  number  of  axial  steps 
Maximum  number  of  axial  steps 
Number  of  implicit  iterations 


Method  for  calculating  effective  radius 

Streamline  ^  Geometry  None 


FIGURE  A-1 1.  SAMPLE  INPUTS  FOR  BOUNDARY 
LAYER  EXECUTION  PARAMETERS 


Enter  the  values  shown  in  the  figure  below  and  click  on  'OK' 
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Execution  Parameters  Modeling  Parameters  |  Output  Parameters  | 

Transition  Parameters — - - — * 

^  Laminar 


oi 


Axial  tronsHion  location 
('  Gradual  tronsHion  factor 

Modeling  Factors - ^ - 


Minimum  pressure  gradient  to  initiate 
shock  seooration 

Streamline  divergence  factor 
Streamline  convergence  factor 
Induced  pressure  limHer 


Momentum  and  energy  thickness  limiter  1 0 


OK 

Cancel 

Apply 

Help 

FIGURE  A- 12.  SAMPLE  INPUTS  FOR  BOUNDARY 
LAYER  MODELING  PARAMETERS 


A.3.5  OUTPUT  CONTROLS 

This  section  is  used  to  specify  what  type  of  output  is  generated,  as  well  as  the  frequency 
of  the  output  (with  respect  to  axial  integration  steps).  Each  of  the  parameters  shown  below  are 
described  in  detail  in  the  Output  Controls  section. 

•  Click  on  'Inputs'  and  then  select  'Output  Control' 

•  Enter  all  of  the  values  shown  in  the  figure  below  and  click  'OK'. 
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f  Surface  Pressures  -  Forces  -  Moments - - - ■ "  “ 

I  C  Print  cell  coordinates 
]  Do  not  print  cell  coordinates 

I  Write  surface  properties  at  every  V  axial  station  1^^  ~ 

j  Write  surface  properties  if  stepsize  is  greater  than  DE1_2A  5000 


HilBlI 


Contour  Data  - - - - - . . . .  — _____ 

Write  contour  properties  at  every  'k'  axial  station  1^00  _ 

and  _ 

atthe  following  user  specified  axial  locations  |0.20  0.22  0.23.50 .51 .52  _ 

(Must  be  in  ascending  order)  _ _ ^ 

Wiite  contour  properties  if  axial  stepsize  is  greater  then  DELZC.  |5000 


-ASCII  Data - - - : - - — - - - 

Output  residual  every  ISKIP  iterations 

Print  crossflow  plane  if  step  number  is  evenly  divisible  by  IPRINT  jsoOO 
Print  N  planes  which  are  evenly  divisible  by  NSKIP  |l 

Print  M  planes  which  are  evenly  divisible  by  MSKIP  1 


FIGURE  A- 13.  SAMPLE  INPUTS  FOR  OUTPUT  CONTROL  SECTION 


A.4  EXECUTION 


RUNNING  A  CASE 

This  section  is  used  to  select  which  cases  are  evaluated  as  well  as  what  type  of  solution 
is  performed. 

•  Click  'Execution'. 

•  Click  on  the  desired  run. 

Note:  Multiple  runs  may  be  selected  in  a  number  of  ways.  The  first  is  to  hold  down 
the  'Ctrl'  key  and  click  on  each  of  the  desired  runs.  Alternatively,  you  can  click 
on  the  first  desired  run,  hold  down  the  ’shift'  key,  and  then  click  on  the  last  run 
to  select  all  cases  in  between.  Finally,  you  can  simply  click  on  the  ’select  AH' 
button  to  run  all  cases. 
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•  Select  the  type  of  run  desired. 

Note;  If  'Full  Run’  is  selected  the  code  will  execute  the  Inviscid  Euler  solver,  write  out 
the  flowfield  at  the  specified  axial  intervals,  and  then  execute  the  boundary  layer 
solver  to  get  an  approximate  viscous  solution.  If  Inviscid  Run'  is  selected,  the 
code  will  execute  the  Inviscid  Euler  solver  and  then  terminate.  No  viscous 
approximation  will  be  calculated.  If  a  viscous  calculation  is  desired  at  later  time, 
simply  click  on  the  'Viscous  Run'  button,  and  select  the  desired  case.  Only  those 
runs  that  have  been  solved  with  the  inviscid  solver  will  be  displayed  when  the 
'Viscous  Run'  button  is  selected.  This  is  because  the  viscous  approximation 
requires  the  inviscid  flowfield  as  a  boundary  condition. 

•  Specify  whether  a  symmetry  plane  exists. 

Note:  If  you  know  a  priori  that  the  geometry  is  non-symmetric  about  the  pitch  plane, 
you  must  select  'Full  Grid  (All  Runs)'.  However,  if  the  geometry  is  symmetric 
about  the  pitch  plane  choose  either  ’symmetry  Plane  (All  Runs)'  if  you  will  not 
be  running  any  asymmetric  cases,  or  alternatively,  choose  ’symmetry  Plane 
(Zero  Roll  and  Yaw)’  to  run  a  symmetry  plane  only  for  those  cases  with  zero  yaw 
or  roll  angles.  Running  with  a  symmetry  plane  will  cut  the  computational  time 
approximately  in  half. 

•  Click  on  'OK'  to  perform  the  desired  calculations. 

•  After  the  runs  have  completed,  you  can  go  back  to  the  'Execution'  section  to  check  on 

the  termination  status  for  each  run.  If  the  'Completed'  column  contains  one 
then  the  inviscid  solution  has  completed.  If  it  contains  '**'  then  both  the  inviscid 
and  viscous  solutions  have  been  completed.  In  addition,  the  'Exit  Status'  column 
should  contain  'Normal  Exit'  if  no  problems  were  encountered.  If  the  code  did 
not  exit  normally  then  an  exit  code  will  be  provided.  The  exit  codes,  along  with 
an  explanation  of  each,  can  be  found  in  the  Exit  Codes  section. 
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HlilEa 


Select  Runs  (Press  Ctrl  key  for  muttipie  selections) - 

Attack  Yaw  Roll 


Run  Exit  Status 


RunOOOi  i 

4000 

0.000 

0.000 

0.000 

Normal  Exit 

Run0002 

4.250 

0.000 

0.000 

0.000 

Normal  Exit 

Run0003 

4500 

0.000 

0.000 

0.000  •: 

Normal  Exit  i 

RunOOO^ 

4.000 

10.000 

0.000 

0.000  - 

* 

Normal  Exit 

A. 

RunOOOB 

4.250 

10.000 

0.000 

0.000  ; 

Jt 

Normal  Exit  ‘ 

RunOOOB 

4.500 

10.000 

0.000 

0.000 

A 

Normal  Exit  , 

Run0007 

4000 

20.000 

0.000 

0.000  ; 

1 

RunOOOS 

4.250 

20.000 

0.000 

0.000  ^ 

r 

Run0009 

4.500 

20.000 

0.000 

0.000 

u 

RunOOl  0 

4000 

30.000 

0.000 

0.000 

sic 

51 05  ; 

f 

RunOOII 

4.250 

30.000 

0.000 

0.000  i 

I 

RunOOl  2 

4.500 

30.000 

0.000 

0.000  : 

Select  AJl 


-  Euler/Boundary  Layer  Run- - — 

r  Full  Run 
Invisdd  Run 

r  Viscous  Run  (Boundaiy  Layer) 


r 


Full  grid/Symmetry  plane - — 

r  Full  Grid  (AJl  runs) 
r  Symmetry  Rane  (AJl  runs) 

Symmetry  Plane  (Zero  yaw  and  roll) 


OK 


. , _ 


Cancel 


in 


FIGURE  A-14.  SAMPLE  SCREEN  SHOT  OF  EXECUTION  CONTROL  PARAMETERS 


A.5  POST-PROCESSING 


OVERVIEW 

The  post-processing  section  of  the  tutorial  is  divided  into  three  sections.  To  execute  the 
plotting  functions  of  the  post-processor,  a  current  version  of  the  Tecplot  plotting  software  must 
be  installed  on  the  user’s  machine.  If  this  software  is  not  installed,  the  user  may  still  examine 
the  ASCII  output  data  generated  by  Zeus**. 
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ASCII  Data  Examine  an  ASCII  data  file  of  the  forces/moments  acting  on 

the  body  as  well  as  the  centers  of  pressure. 

Single  Run  Plots  Line  plots  of  the  surface  pressures,  forces,  and  moments. 

Vector  and  streamline  plots.  Contour  plots  of  a  variety  of 
flowfield  parameters  (Mach  number,  density,  temperature, 
...etc). 

Multi-Run  Plots  Comparative  line  plots  of  forces/moments  versus  a  variety  of 

independent  variables  (Mach  number,  angle  of  attack,  yaw 
angle,  and  roll  orientation). 

A.5.1  ASCn  FORCES  &  MOMENTS 

This  section  is  used  to  examine  the  ASCII  output  data  from  the  Zeus**  runs.  The  data 
includes  the  forces,  moments,  and  centers  of  pressure  for  both  the  entire  geometry  as  well  as 
each  individual  zone  edge.  An  in-depth  explanation  of  the  ASCII  data  can  be  found  in  the 
ASCII  Forces  &  Moments  section  of  the  users'  manual. 

•  Click  on  'Post-Processor'  and  then  select  'ASCII  Forces  &  Moments' 

•  Select  'Totals'  to  display  the  total  forces/moments  acting  over  the  entire  geometry. 
Alternatively,  'Edge  Values'  could  be  selected  to  display  the  individual  components  of 
the  forces/moments  that  constitute  the  total. 

•  Microsoft  WordPad  will  open  and  display  the  ASCII  (forcesmoments.txt)  version  of 
the  output  data.  The  data  can  either  be  printed  or  saved  to  another  file  using  this 
program.  When  finished  examining  the  data,  close  the  WordPad  program  to  return  to 
Zeus**. 


A.5.2  SINGLE-RUN  PLOTS 

This  section  is  used  to  generate  graphs  (both  line  and  contour)  of  the  forces,  moments, 
centers  of  pressure,  and  a  variety  of  flowfield  parameters  (density,  Mach  number,  pressure, 
etc...).  Comparison  plots  of  the  various  runs  are  generated  in  the  Multi-Runs  section. 

■  Click  on  'Post-Processor'  and  then  select  ’single  Run  Plots'. 

■  Sel^t  'Plot  Surface  Pressure'.  A  window  will  then  be  displayed  which  lists  all 
completed  runs.  Select  the  desired  run  and  click  on  'OK'  to  generate  the  graph.  A 
detailed  description  of  the  data  in  the  surface  pressure  plot  can  be  found  in  the  users' 
manual  in  the  Surface  Pressure  section. 

■  Exit  the  plotting  package  to  return  to  Zeus**. 

■  Click  on  'Post-Processor'  and  then  select  ’single  Run  Plots'. 

■  Select  'Plot  Forces  &  Moments'  and  then  choose  either  'Totals'  or  'Edge  Values'  to 
generate  the  corresponding  graph.  After  the  completed  runs  are  displayed,  select  the 
desired  run  and  click  on  'OK'  to  generate  the  graph.  A  description  of  the  data  in  the 
graph  is  provided  in  the  users'  manual  in  the  Forces  &  section. 
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■  Exit  the  plotting  package  to  return  to  Zeus**. 

■  Click  on  'Post-Processor'  and  then  select  ’single  Run  Plots'. 

-Select  'Contour  Plots’  and  then  'Choose  Plot  Attributes’.  Enter  the  desired 
dimensionalization  parameters  for  the  contour  plots,  select  the  desired  variables,  and 
click  on  'OK'.  The  above  parameters  are  explained  in  detail  in  the  ContourA^ector 
section  of  the  user’s  manual. 

■  Return  to  the  'Contour  Plots'  menu,  however,  this  time  click  on  Plot  Data  to  generate 
the  contour  graph.  After  the  completed  runs  are  displayed,  select  the  desired  run  and 
click  on  'OK'  to  generate  the  graph. 

■  Exit  the  plotting  package  to  return  to  Zeus**. 


A.5.3  MULTI-RUN  PLOTS 

This  section  is  used  to  make  comparison  plots  of  different  runs.  The  options  include 
comparisons  of  all  forces,  moments,  and  centers  of  pressure  with  a  choice  of  dependant 
variables  (Mach  number,  angle  of  attack/yaw/roll). 

■  Click  on  'Post-Processor'  and  then  select  'Multi-Run  Plots'. 

■  Select’Choose  Plot  Attributes'  and  specify  the  values  to  be  plotted  as  well  as  the 
independent  variable.  Re-fer  to  the  Multi-Plot  run  section  of  the  user’s  manual. 

■  Return  to  the  'Multi-Run  Plots'  menu,  however,  this  time  click  on  Tlot  Data’  to 
generate  the  comparison  plot  for  all  executed  runs.  If,  during  the  tutorial,  only  a 
single  run  was  executed,  then  the  line  plot  displayed  will  be  a  point.  In  order  to  make 
comparative  plots  of  different  runs,  return  to  the  execution  section,  and  run  additional 

cases. 

■  Exit  the  plotting  package  to  return  to  Zeus**. 
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